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ABSTRACT 

Exact  high  temperature  series  expansions  have  been 
obtained  for  various  thermodynamic  quantities  for  the 
spin  one  half  XY  model  of  ferromagnetism  or  of  a  quantum 
lattice  fluid. 

The  static  quantities  calculated  include  the  initial 
parallel  susceptibility  and  the  fourth  order  fluctuation 
both  for  the  face  centered  cubic  and  for  the  triangular 
lattice.  Analysis  of  the  series  for  these  two  quantities 
by  standard  ratio  and  Fade  approximant  methods  yields 
estimates  for  the  critical  temperature  variable,  Kc=J/kTc' 
and  the  critical  indices  y  and  y2  respectively.  For  the 
face  centered  cubic  lattice  Kc=  0.221±0.001,  y=1.32+0.04 
and  y 2=  4.6410.10.  For  the  triangular  lattice  Kc=° • 67±0 . 01 , 
y=l .5010.02  and  y2=  5.410.5.  These  results  coupled  with 
scaling  theory  yield  all  critical  indices  for  the  spin  one 
half  XY  model.  On  the  triangular  lattice  the  results 
obtain  if  we  assume  a  phase  transition  but  do  not  consti¬ 
tute  conclusive  evidence  of  the  existence  of  such  a 
transition.  The  results  for  the  face  centered  cubic 
lattice  compare  very  well  with  experimental  results  for 

the  X  transition  in  liquid  helium. 

The  dynamical  calculations  have  been  done  for  two 
quantities,  the  autocorrelation  and  the  frequency 


' 


dependent  parallel  susceptibility,  both  on  the  face 
centered  cubic  lattice.  The  autocorrelation  can  be 


in 


expanded  in  a  power  series  in  time.  We  have  calculated 

the  first,  second,  third,  fourth,  sixth  and  seventh 

coefficients  as  power  series  in  T  ^ .  We  find  that  the 

first,  second  and  third  coefficients  diverge  at  K  with 

c 

exponents  y^  =  0.1010.01,  y^  =  0.24±0.05  and  y^  =  0.3510.05 
respectively . 

The  frequency  dependent  susceptibility  can  be  expanded 

in  a  power  series  in  the  reciprocal  of  the  frequency,  m  ± . 

Coefficients  of  odd  powers  vanish;  coefficients  of  even 

powers  have  been  obtained  as  power  series  in  the  reciprocal 

temperature  variable,  K.  We  find  that  the  second  moment 

diverges  at  K_  with  exponent  yn=  0.10+0.01  and  the  fourth 
— — — — — —  c  ** 

moment  vanishes  with  exponent  y,=  -0.86+0.05.  Assuming 

_  y  As 

that  a  relation  x(a)/T)  =  [(TC“T)/TCJ  f  (to/  [  (T-Tc)  /Tc]  ) 

holds  near  T  for  all  w  we  find  that  A  =  0.5810.10.  It 

c 

is  important  to  observe  that  Ag  <  y • 
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CHAPTER  I 


INTRODUCTION 

This  thesis  describes  a  study  of  the  thermodynamical 
properties  above  the  critical  temperature  for  the  spin  h 
XY  model.  The  method  of  exact  high  temperature  series 
expansion  is  used.  The  XY  model  is  quantum  mechanical  and 
yet  calculations  with  it  are  relatively  easy.  The  spon¬ 
taneous  magnetization,  which  is  the  order  parameter,  lies 
in  the  xy  plane.  Therefore  the  order  parameter  in  the  XY 
model  does  not  commute  with  the  Hamiltonian  and  can 
exhibit  dynamical  behaviour.  This  is  not  the  case  in 
the  Ising  and  Heisenberg  models,  where  the  order  parameter 
is  constant  in  time  and  changes  in  the  Hamiltonian  have 
to  be  introduced  to  allow  dynamical  behaviour  of  the  mag¬ 
netization.  Accordingly  this  thesis  is  divided  into  two 
parts .  The  first  deals  with  static  critical  phenomena  and 
the  second  treats  dynamical  critical  phenomena. 

In  Section  2.2  the  XY  model  is  introduced.  We  shall 
briefly  summarize  here  the  reasons  for  the  interest  in 
this  model.  In  1956  Matsuda  and  Matsubara  proposed  the 
anisotropic  Heisenberg  model  as  a  possible  model  for  an 
interacting  Bose  fluid.  The  XY  model  represents  a  Bose 
fluid  with  hard  core  repulsive  interaction  only. 
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The  classical  anisotropic  Heisenberg  model  has 
been  studied  by  Jasnow  and  Wortis  (1968)  for  various 
values  of  the  anisotropy  parameter.  They  concluded 
that  the  critical  behaviour  for  the  classical  model 
depends  only  on  the  dimensionality  of  the  lattice  and 
the  symmetry  of  the  ground  state.  That  is,  the  critical 
indices  are  constant  over  ranges  of  the  anisotropy 
parameter  and  change  discontinuously  when  the  symmetry 
of  the  ground  state  changes.  The  XY  behaviour  persisted 
as  long  as  the  ground  state  was  invariant  under  rotation 
in  the  plane.  This  is  true  when  the  XY  part  of  the  inter¬ 
action  is  stronger  than  the  Z  or  Ising  part. 

Jasnow  and  Wortis  studied  only  the  spin  infinity, 
or  classical,  case  but  the  conjectured  behaviour  is 
supported  for  the  quantum  mechanical  case  by  the  work 
of  Barouch  and  McCoy  (1969) .  Barouch  and  McCoy  found, 
for  the  anisotropic  XY  model  in  one  dimension,  that  the 
nature  of  the  decay  of  the  spin  correlations  changed 
when  changes  in  the  anisotropy  parameter  and  in  the 
applied  field  affected  the  symmetry  of  the  ground  state. 

It  seems  therefore  reasonable  to  expect  that  the  critical 
indices  of  the  spin  ^  anisotropic  Heisenberg  model  will 
be  equal  to  the  indices  of  the  XY  model  as  long  as  the 


■ 
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- 
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XY  interaction  is  dominant.  The  indices  we  obtain  are 
expected  then  to  fit  the  A  transition  of  liquid  helium, 
and  some  insulating  magnetic  crystals  for  which  >>  J, (  , 

that  is  the  transverse  interaction  is  much  stronger  than 
the  parallel. 

The  quantities  usually  expanded  in  high  tempera¬ 
ture  series  are  free  energy,  the  susceptibility  and  the 
even  derivatives  of  the  free  energy  with  respect  to  the 
field  at  zero  field.  Betts,  Elliott  and  Lee  (1970) 
calculated  the  free  energy  and  fluctuation  series  for 

a  variety  of  lattices.  Because  the  magnetization  and 

1  2 

the  Hamiltonian  do  not  commute,  the  fluctuation  ,  ^-<Mx>, 
is  not  equal  to  the  susceptibility.  Falk  and  Bruch  (1969) 
proved  however  that  the  critical  index  of  the  fluctuation 
is  equal  to  y,  the  index  of  the  susceptibility.  The 
proof  involved  only  rigorous  thermodynamical  inequalities. 
The  equality  is  only  asymptotically  true  close  to  the 
critical  temperature.  It  is  therefore  of  interest  to 
calculate  the  susceptibility  to  see  numerically  how  the 
coefficients  approach  those  of  the  fluctuation  series. 

From  the  fluctuation  series  on  the  triangular  lattice 
no  conclusions  could  be  drawn  as  to  the  existence  of  a 


■ 
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transition  and  it  was  hoped  that  the  static  suscepti¬ 
bility  would  behave  better. 

The  other  static  quantity  calculated  was  the  fourth 
order  parallel  fluctuation 


„  1  r  3  ,2  2  1  ..4  , 

Y2  =  N  {  2  "“x"  "  2  <Mx"  } 


By  similar  arguments  to  those  used  for  the  fluctuation 
by  Falk  and  Bruch  (1969) ,  diverges  as  the  fourth 
derivative  of  the  free  energy  with  respect  to  the  field 
in  the  x  direction,  thereby  giving  us  an  estimate  for  the 
gap  parameter.  Now  two  critical  exponents  and  scaling 
yield  estimates  for  all  critical  indices. 

Before  leaving  the  static  quantities  we  should 
mention  that  the  limit  at  zero  frequency  of  the  dynamical 
susceptibility  is  the  quasistatic  (called  also  adiabatic 
after  Kubo)  susceptibility  while  ours  is  the  isothermal. 
They  are  equal  when  the  system  is  ergodic.  This  question 
arose  in  the  dynamical  calculations  for  the  Ising  model 
of  Allan  and  Betts  (1968)  and  Essam  and  Garelick  (1968) . 
We  shall  discuss  ergodicity  in  Chapter  V  but  mention 
here  the  conclusion  that  in  the  XY  3-dimensional  case 
we  believe  the  two  are  equal. 


, 


, 
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Dynamic  critical  phenomena  are  at  a  much  earlier 
stage  of  theoretical  investigation  than  static  critical 
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phenomena . 

Scaling  laws  have  less  success  with  static  corre¬ 
lations  than  they  had  with  the  thermodynamical  quantities. 
For  the  dynamical  correlations  the  assumptions  used  for 
the  static  correlations  are  kept  and  generalized,  there¬ 
fore  scaling  is  on  a  weaker  basis. 

There  is  no  exactly  soluble  model  with  physical 
critical  behaviour  that  can  be  used  to  check  dynamical 
theories  the  way  the  2d  ising  model  is  used  in  static 
theory.  The  only  soluble  model  is  the  one  dimensional 
XY  model.  Niemeijer  (1967)  ,  Barouch,  McCoy  and  Dresden 
(1969)  and  Suzuki  (1969)  have  calculated  the  time  de¬ 
pendent  correlations  and  magnetization  of  the  one 
dimensional  XY  model .  This  model  has  the  disadvantage 
of  being  non  ergodic .  The  limit  t  -»•  °°  yields  a  magneti¬ 
zation  which  is  not  the  equilibrium  magnetization  in  the 
2  direction.  The  magnetization  in  the  x  direction  can 
be  expected  to  be  ergodic  and  of  interest  but  cannot  be 
calculated  exactly  with  present  methods. 

The  kinetic  Ising  model  in  2  dimensions  was  inves¬ 
tigated  by  Suzuki  (1969)  using  high  temperature  series 
expansions.  Suzuki  found  that  the  critical  index  of 
slowing  down  is  different  from  the  index  of  the  static 


. 

- 

- 
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susceptibility.  Dynamical  molecular  field  theory  (Van 
Hove  (1959),  Suzuki  and  Kubo  (1968))  predicts  equality 
for  the  two  indices. 

Recently  McFadden  and  Tahir  Kheli  (1970)  calculated 
the  zeroth,  second  and  fourth  moments  of  the  spectral 
function  as  functions  of  the  momentum  transfer.  Only 
the  leading  term  in  T  ^  was  obtained. 

We  have  calculated  the  frequency  dependent  suscep¬ 
tibility  x (0/  •  We  have  obtained  a  double  series  in 

-1  -1 

powers  of  T  and  oj  .  Analysis  of  the  series  is  hampered 
by  the  fact  that  coefficients  of  powers  of  T  ^  are  infinite 
series  in  powers  of  w  1  and  vice  versa.  We  find  that  the 
second  moment  of  the  susceptibility  diverges  at  T  >  Tc . 

For  the  kinetic  Ising  model  and  the  isotropic  Heisenberg 
model  all  moments  are  finite,  as  was  rigorously  proved  by 
Suzuki  (1969)  and  Mermin  and  Wagner  (1966)  respectively. 
Though  the  proofs  could  not  be  generalized  to  the  XY  model 
the  divergence  was  not  expected . Assuming  scaling  we 
estimate  the  critical  index  of  slowing  down, and  we  find 
that  in  this  case  it  is  speeding  up.  We  calculated  the  auto¬ 
correlation  and  with  a  scaling  assumption  estimate  its 


relaxation  time. 


' 
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CHAPTER  II 

STATIC  CRITICAL  PHENOMENA  AND  THE  XY  MODEL 

2 . 1  Critical  phenomena  and  critical  exponents 

Phase  transitions  have  been  studied  by  physicists 
since  the  day  man  discovered  that  water  boiled  and  froze. 

A  great  variety  of  such  phenomena,  while  of  completely 
different  origin,  have  a  very  similar  quantitative  behaviour 
and  can  be  studied  using  the  same  mathematical  models. 

Some  examples  of  critical  phenomena  are  those  which 
occur  in  monomolecular  fluids  near  the  critical  point,  in 
liquid  helium  near  the  A  temperature,  in  binary  alloys  near 
the  order-disorder  transition  and  in  ferromagnets  in  the 
vicinity  of  the  Curie  temperature. 

In  all  these  systems  by  varying  a  thermodynamical 
parameter  we  reach  a  point  where  there  is  no  phase  transi¬ 
tion,  only  one  homogeneous  phase  exists  and  changes  such  as 
in  density  or  magnetization  are  continuous. 

For  all  those  phenomena  a  long  range  order  parameter 
R  can  be  defined.  At  T=0  there  is  absolute  order  in  the 
system.  All  spins  are  lined  up  parallel  to  each  other  in 
a  ferromagnet,  all  He^  atoms  are  in  the  ground  superfluid 
state  etc.,  so  R=Rq;  at  infinite  temperature  R=0  as  complete 
disorder  must  occur.  In  the  systems  with  which  we  are 
concerned  however  the  long  range  order  disappears  suddenly 
at  and  above  some  temperature  defined  as  the  critical  point. 


. 
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whereas  for  systems  which  do  not  exhibit  phase  transitions 
R  decays  to  zero  smoothly  and  only  at  T=°° ;  for  the  transi¬ 
tions  with  which  we  deal  the  behaviour  of  R  is 


For  example,  the  order  parameter  for  a  ferromagnet  is  the 
magnetization,  and  for  superfluid  it  is  the  density  of  the 
superfluid  fraction. 

In  the  vicinity  of  this  critical  point  there  are 
large  fluctuations  of  all  thermodynamic  quantities  which 
can  become  macroscopic  and  the  relevant  derivatives  of  the 
thermodynamic  potential,  such  as  specific  heat  and  suscep¬ 
tibility,  diverge .  Such  macroscopic  fluctuations  cannot 
probe  deeply  into  the  microscopic  details  of  the  forces 
involved  and  this  is  the  reason  of  the  quantitative  simi¬ 
larity  among  such  diverse  phenomena. 

The  determination  of  the  asymptotic  laws  which 
describe  the  approach  to  the  critical  point  has  been  the 
main  problem  in  the  study  of  critical  phenomena.  It  may 
be  that  only  the  dimensionality  and  symmetry  of  the  system 
are  needed  to  predict  the  asymptotic  laws  of  the  appro¬ 
priate  variable  .  This  is  a  unifying  pleasing  hypothesis 
(Jasnow,  Wortis  1968;  Kadanoff  1970). 
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(a)  Magnetic  systems 

Ferromagnets  are  defined  as  those  crystals  which 
have  a  spontaneous  magnetization  below  the  Curie  tem¬ 
perature  T  ,  that  is  which  possess  a  magnetic  moment  even 
o 

in  the  absence  of  an  applied  magnetic  field  H. 


M  (T)  =  lim  M  (H ,  T)  (2.1.1) 

°  H->0  + 

Above  is  the  paramagnetic  region  and  the  magnetization 
varies  continuously  with  H  through  the  line  H=0. 

The  magnetic  exponent  $  is  defined  by 

M  (T)  ~  (T  -  T  ) ^  (2.1.2) 

o  c 

as  T  approaches  T_.  This  behaviour  is  well  verified  by 

o 

the  experimental  data  of  Heller  and  Benedek  (1965)  and 
Senturia  and  Benedek  (1966)  on  EuS  and  CrBr^- 
The  isothermal  zero  field  susceptibility 


f— ) 

XT  l3H  H=0 


(2.1.3) 


diverges  as  the  temperature  approaches  T( 

T  >  T 


(T-Tc)  Y 


X 


T 


(2.1.4) 


(T-T  ) 
c 


■y 


T  <  T 


This  relation  defines  the  critical  exponents  y  and  y'  and 
again  experiments  bear  out  this  assumed  behaviour  for  a 


■ 
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variety  of  materials: Ni ,  Fe,  Gd ,  YtFeO^  (see  Miedema  et  al 
(1963) ) .  The  critical  index  6  is  defined  on  the  critical 
magnetic  isotherm 

H  ~  I m | 6  T  =  T  .  (2.1.5) 

1  i  c 

The  specific  heat  critical  exponents  a  and  a'  are  defined 

(2.1.6) 

singularity , for  example 

I 

A  can  be  defined  (for 
n 

(2.1.7) 

where  F  is  the  free  energy,  and  similarly  is  defined 

for  T  >  T  . 

c 

For  n  =  1,  (2.1.7) does  not  define  new  exponents  as 
they  are  related  to  a ,  3 /  Y '  by 

A1  =  2  “  a'  “3 
A^  =  3  +  Y' 

Above  T  odd  derivatives  vanish  for  H=0  because  M ( -H , T) =-M (H , T) 
c 

A  and  A  cannot  be  measured  experimentally  for  n>  2. 
n  n 


oy 


'H=0 


(T-T  ) 
c 


■a 


(T-T  ) 
C 


-a 


T  <  T 


T  >  T 


a=0  is  by  definition  a  logarithmic 
EuO  (Teaney  1966)  . 

One  more  series  of  exponents 


T  <  T  )  by  the  relation: 
c  J 


(— ) 

'  9H'  T 


(T  -T) 
c 


’An  /5n  1F\ 
9Hn_1/  T 


■ 
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The  exponents  v,v !  and  p  refer  to  the  behaviour  of 
the  pair  correlation  function  V (r)  and  correlation  length 
E,  observed  in  neutron  scattering^  first  explained  by  Van 
Hove  (1945) 


€ 


'(T-TJ  V’ 

,<T-Tcrv 


T  <  T 

c 


T  >  T 

c 


(2.1.8) 


r 

c 


(r) 


1 

d-2+g 

r 


(2.1.9) 


where  d  is  the  dimensionality  of  the  system. 


The  spin  correlation  function  is  defined  by 
<ca  q^>  —  <  Sa><  S3> 

r  (r,H,T)  =  °  - - 5 - —  a,B=x,y,z  (2.1.10) 

P  j  S  (S  +  1) 


and  T  is  the  correlation  function  in  the  preferred  direc¬ 
tion  at  the  critical  point.  The  quasi  elastic  scattering 
intensity  for  a  momentum  transfer  q  is: 


I  (q) 
I°(q) 


1  +  E  e1-'-  f (r) 


(2.1.11) 


The  coherent  scattering  at  zero  field  below  T^  is 

I  ,  (0,T)  ~  M2 (T)  ~  (T-T)23  (2.1.12) 

coh  o  c 

which  gives  another  experimental  way  of  measuring  3 • 
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The  zero  angle  scattering  is  proportional  to  the  suscepti¬ 
bility  , 

1(0)  ~  xT  ~  (T-Tc)"Y  H  =  0  (2.1.13) 


T  >  T 

c 

Experiments  were  performed  by  Passel  et  al  (1965)  ,  Bailey 
et  al  (1967)  on  iron  to  measure  y. 

From  (2.1.9  )  and  (2.1.11)  we  see 


I  ...  (k)  ~  -= - 

critical  -  ,  2—7~i 


as  k  ->  0 


(2.1.14) 


Jacrot  et  al  (1962)  ,  Passel  et  al  (1965)  measured  g  this 
way . 

The  spin  correlation  function  decays  as 

r  (r ,  T)  -  sJh  .  (2.1.15) 

- '  r 


The  correlation  length  is  defined  by 


K  =  K 


-1 


(2.1.16) 


-1  2 

k  (T)  measures  the  slope  of  I  (k,T)  against  k  as  k  0 


-2 


fr  T (r)  dr 
1  +  V-1/r  (r)  dr 


(2.1.17) 


Near  T  we  expect  K ^  ~  k  if  near  the  critical  point  only 
one  divergent  temperature  dependent  correlation  length 
exists.  The  quantity  v  that  is  measured  is  the  one  per¬ 
taining  tO  K-^. 
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This  was  a  short  summary  of  the  definitions  of 
critical  indices  for  uniaxial  f erromagnets ;  our  interest 
is  really  in  ferromagnets  with  a  plane  of  easy  magnetiza¬ 
tion  but  applying  a  small  field  in  a  direction  lying  in 
the  plane  enables  us  to  treat  that  direction  as  the 
preferred  axis  of  magnetization  and  all  definitions  still 
apply. 


(b)  Superfluid  helium 


At  first  sight  the  superfluid  transition  seems  rather 
different  from  the  transition  in  a  uniaxial  ferromagnet 
because  of  the  existence  of  the  X  line  (see  Figure  2.2) . 

The  confusion  arises  due  to  regarding  the  pressure  on  the 
superfluid  as  the  analogue  of  the  magnetic  field  and  the 
density  as  the  analogue  of  the  magnetization.  The  correct 
analogy  is  between  the  superfluid  order  parameter  (the 
square  root  of  the  superfluid  density)  and  the  magnetic 
order  parameter  (the  axial  magnetization) .  The  generalized 
force  conjugate  to  the  superfluid  order  parameter  and 
analogous  to  the  axial  magnetic  field  is  unfortunately 
not  physically  realizable.  The  analogues  of  the  density 
and  pressure  are  more  correctly  in  the  uniaxial  magnet, 
the  perpendicular  magnetization  and  the  perpendicular 
magnetic  field.  We  shall  expect  therefore 


no(t)  -  (T-Ta)2B 

n  is  defined  by 
o 


n 

o 


=  lim  <H,i(r)H'(r')> 
(r-r '  ) 


(2.1.18) 


2 


(2.1.19) 


. 


- 


' 
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Y (r)  being  the  wave  function  in  second  quantization 
formulation.  Most  simple  theories  of  superfluidity 
assert 


p  being  the  density  of  the  superfluid.  Measurements 
s 

by  Clow  and  Reppy  (1966)  and  Tyson  and  Douglass  (1966) 

show  p  (T)  ~  (T  -T) C  T< T  with  £  -  0.666  ±  0.1, 

s  X  1 

i 

that  is  3  ~  ^  •  T^e  sPecific  heat  at  constant  pressure 
shows  a  logarithmic  singularity  (Figure  2.3).  Experimental 
results  by  Fairbank  and  Kellers  (1966),  and  more  accurate  results 
by  Ahlers  (1969),  give  a=  0.000  ±  0.003  according  to  one 
interpretation  of  the  data  and  a  =  -0.005  ±  0.005  according 
to  another. 


C  ~  -A  In  (1  -  =-)  +  B  T  >  T, 
P  T,  ' 


(2.1.20) 


Henkel,  Smith  and  Reppy  (1969)  obtained  vh=  0.67  ±  0.04 
by  measuring  the  superfluid  healing  length  of  thin  films 
which  by  the  assumption  of  one  relevant  coherence  length 

near  T ^  is  the  same  as  v. 

Once  the  Hamiltonian  is  given  the  problem  is  only 

to  evaluate  the  partition  function 

-3# 


Z(T,N,V)  =  Tr  e 


N 


(2.1.21) 


. 


, 


' 


and  then 
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F  (T ,  V) 
kT 


lim 

V-*-00 


1 

N 


In 


Z 


(2.1.22) 


and  all  thermodynamical  quantities  follow  by  the  usual 
relations.  The  problem  of  going  to  the  thermodynamic 
limit  in  order  to  have  any  critical  behaviour  is  at  least 
partially  solved  by  the  rigorous  proofs  of  Ruelle  (1963) 
and  others  who  showed  that  thermodynamical  properties 
are  the  same  whether  obtained  by  canonical  or  grand 
canonical  ensemble  and  do  not  depend  on  the  shape  of  the 
sequence  of  volumes  one  uses.  Still  the  order  in  which 
limits  are  taken  can  make  a  difference  and  caution  is 
required . 


. 
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Figure  2 . 1 

Susceptibility  in  a  Ferromagnet 


Figure  2 . 2 

4 

Schematic  Phase  diagram  of  He 
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Figure  2.3 

Specific  heat  of  liquid  helium 


2 . 2  The  XY  model  and  its  relation  to  magnetic  systems 
and  quantum  fluids 


18 


The  first  attempt  to  build  a  model  for  ferromagnetism 
was  the  Weiss  molecular  field  theory  (1907) .  The  model 
assumes  that  each  particle  in  the  system  interacts  equally 
with  all  other  particles  in  the  system.  In  spite  of  this 
rather  unrealistic  assumption  the  results  were  in  some 
agreement  with  experiment  at  least  in  predicting  a  critical 
point,  a  diverging  susceptibility  (with  6=1)  and  a  dis¬ 
continuous  specific  heat. 

The  first  realistic  model  to  be  treated  was  the 
Ising  model  with  first  neighbour  interactions.  The  spins 
are  fixed  on  a  lattice  and  the  Hamiltonian  describing  the 
interaction  is 

dt  =  -  J  E  qZ  aZ-mHZ  Z  qZ  (2.2.1) 

..in  .1 

<il>  J  l 

where  <ij>  denotes  sum  over  pairs  of  nearest  neighbours 
cZ  is  the  Pauli  spin  matrix  and  m  is  the  magnetic  moment 
of  the  atoms . 

It  was  solved  for  the  spin  ^  case  for  a  one  dimen¬ 
sional  chain  by  Ising  (1920),  and  by  Onsager  (1944)  for  a 
square  lattice.  Series  expansions  gave  good  estimates  for 
the  critical  indices  in  3  dimensional  lattices.  Magnetic 
systems  were  found  to  agree  with  the  predictions  of  this 
model,  namely  the  cobalt  tutton  salts. 


The  8  brass 
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binary  alloy  seems  to  be  a  good  example  for  this  model 
(Als-Nielsen  1969)  . 

The  Heisenberg  model  is  based  on  the  exchange 
theory  of  ferromagnetism.  The  Heisenberg  Hamiltonian  is 

=  -2 J  £  Si.  St  -m  HZ£  S*  •  (2.2.2) 

<  i  j  >  J  i 

Materials  like  EuO  and  EuS  satisfy  both  assumptions  of 
the  model ,  localization  of  spins  and  isotropy  of  inter 
action,  and  series  expansion  methods  gave  estimates  for 
the  critical  indices  (Rushbrooke,  Wood  (1958) ,  Domb,  Wood 
(1964),  Baker,  Gilbert,  Eve  and  Rushbrooke  (1967)). 

The  natural  extension  of  the  Heisenberg  model  is 
a  model  in  which  the  assumption  of  isotropy  is  dropped. 
Mathematically  this  causes  great  complications,  usually 
the  most  general  model  treated  is 


Ji  = 


s  {J ,  (s*s*  +  sM)  +  j 

<ij> 


SZSZ}-  m  H 
i  J 


£ 

i 


i 


(2.2.3) 


which  is  called  the  anisotropic  Heisenberg  model  even 
though  the  anisotropy  is  not  complete.  The  Ismg  modex 

is  the  extreme  case  of  =  0*  The  XY  model  is  the  other 
extreme  J  =0.  The  XY  hamiltonian  is 


xy 


-2  £  J(S^  +  S^) 

<ij> 


(2.2.4) 


( 
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first  introduced  by  Lieb,Schultz  and  Mattis  (1961). 

The  work  done  on  the  anisotropic  Heisenberg  model  was 
to  calculate  yZZ which  for  the  limiting  case  of  the  XY 
interaction  is  not  expected  to  diverge,  whereas  this 
thesis  deals  with  a  field  applied  in  a  direction  in  the 
plane  of  interaction. 

A  strong  inducement  to  study  this  model  was  the 

work  of  Jasnow  and  Wortis  on  the  classical  anisotropic 

Heisenberg  hamiltonian.  The  classical  limit  is  obtained 

-*■  s 

by  taking  the  limit  S  00  but  -TqT  ^  f  •  This  gives  an 
interaction  of  commuting  spins  which  behave  like  classical 
vectors.  Jasnow  and  Wortis  studied  this  Hamiltonian  for 
various  anisotropies  on  different  lattices  by  series 
expansions.  Their  results  suggest,  as  mentioned  in  2.1, 
that  only  the  symmetry  of  the  ground  state  and  the  dimen¬ 
sionality  of  the  lattice  affect  the  critical  exponents. 

In  particular  as  long  as  the  interaction  has  the  symmetry 
of  the  XY  hamiltonian,  that  is  as  long  as  |jg_l>|j]|  i  /  f 
critical  behaviour  of  the  system  is  that  of  pure  XY 
interaction.  The  indices  change  only  when  there  is  a 
change  in  the  symmetry / and  then  they  change  discontinuously . 
The  same  result  is  indicated  by  Barouch  and  McCoy  (1970) 
for  the  one  dimensional  anisotropic  XY  model  (see  2.3). 

While  the  pure  XY  model  is  mathematically  convenient,  it 
has  no  possible  hope  of  describing  a  physical  system.  However 


' 

■ 
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an  anisotropic  Heisenberg  model  with  the  XY  symmetry  in 
which  Jj~  >>  Jjj  is  quite  possible. 

It  was  suggested  by  Huiskamp  (1967)  that  materials 
with  strongly  anisotropic  crystalline  fields  like 
Gd  2  (SO^)  2*  wouid  be  likely  to  have  an  XY  like 

behaviour.  It  was  shown  by  Betts,  Elliott  and  Lee  (1970) 
that  an  insulating  crystal  composed  of  high  half  odd 
spins  with  a  strong  axial  crystalline  field  behaves  much 
like  a  spin  %  XY  system.  The  interaction  being 


$1  =  +  D  E  (S^)2 
i 


J 


E 

<ij> 


-* 


S  .  .  S  . 
i  : 


d  >>  j  >  o 
(2.2.5) 


2 

the  strong  crystalline  field  will  cause  the  S  '  =  ±h  level 

to  be  highly  populated  in  the  critical  region.  Looking 

now  only  at  the  populated  subspace  S  =  ±h  the  ratio  of 

the  SXSX  +  sYsY  part  of  the  interaction  compared  to  the 
l  j  l  j 

2 

S.S.  part  will  be  about  S'  to  1. 
i  j  ^ 

In  such  a  material  we  have  a  dominant  transverse 
coupling,  the  spontaneous  magnetization  vector  M  m  zero 
field  will  lie  in  the  xy  plane.  Introduction  of  a  field 
in  the  z  direction  will  tend  to  rotate  the  magnetization 
towards  the  z  axis  for  large  enough  fields  and  we  end 
with  a  paramagnetic  phase  in  the  xy  plane.  The  expected 


' 


22 


phase  diagram  (Fisher  1968)  is 


where  the  broken  line  indicates  the  transition  may  not 
be  second  order  any  more/ and  we  have  a  A- like  transition 
line.  Introducing  a  field  in  the  x  direction  will  cause 
the  magnetization  to  be  in  the  x  direction  even  in  the 


limit  of  H  0. 

H,io° 


H.--0 


nx 


TT  T 

—  > 

X  X  'i' 

T 

Our  calculation  is  a  series  expansion  for  T  >  0. 

We  approach  Tc  from  above  on  the  Hx=0  line.  Using  the 
equalities  predicted  by  scaling,  we  can  obtain  the 
behaviour  on  other  lines,  for  example  the  approach  to 

H  =  0  on  a  T=T  line, 
x  c 

Matsuda  and  Matsubara  (1956)  introduced  the  aniso¬ 
tropic  Heisenberg  model  as  a  possible  model  for  a  quantum 
lattice  fluid.  In  the  following  pages  we  bring  Matsuaa  and 
Matsubara ' s  derivation.  The  model  for  a  classical  fluid  as 


' 
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proposed  by  Yang  anu  Lee (1952)  assumes  that  each  atom  occupies 
a  cell  i  with  a  probability  t^  which  is  either  1  or 
0.  Multiple  occupancy  is  thus  excluded  and  the  grand 
partition  function  for  the  system  obtained  is: 


E  «  £  exp  £  t .  -  £  f3cj)  (r  .  -r  .  )  t .  t  .]  (2.2.6) 

t=0 , 1  i  1  ij  1  9  1  D 


li  being  the  chemical  potential  and  (^(r^-ru)  the  potential 
between  atoms  in  cells  i  and  j .  Taking 


4>(r±-r.)  = 


' AJ  for  i  and  j  nearest  neighbours 

(2.2.7) 

0  otherwise 


the  model  reduces  to  the  classical  Ising  model. 

For  a  quantum  fluid  we  work  in  the  framework  of 
second  quantization  and  define  creation  and  annihilation 

JL 

operators  a^  and  a_^  respectively  on  each  lattice  site  i. 
Since  we  are  dealing  with  a  Bose  gas 


=  0 


(2.2.8) 


for 

i=j 


i^j  • 


In  order  to  exclude  multiple  occupation  for 


(a^ ,  at}  =  a^a^  +  ~  1 

(at,  at}  =  (a.,  a.}  =  0 

i'  j  ii 


(2.2.10) 


, 
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and  this  ensures  that  the  number  operators 

N .  =  a  .  a  . 

1  11 


will  have  eigenvalues  1,  0  only.  These  commutation 
relations  are  satisfied  by 


a 


t 

i 


+  i 

l 


x  .  y 
a  .  =  a  .  -  i  a  . 

ill 


(2.2.11) 


The  following  potential  energy  <j>  is  considered: 


for  two  atoms  on  same  lattice  point 
(forbidden  by  (2.2.10)) 
for  nearest  neighbours 
otherwise  . 


Therefore 


< P 


(2.2.12) 


As  to  kinetic  energy  it  is  assumed  that  atoms  can  move 

only  to  vacant  nearest  neighbour  sites.  Such  a  transi- 

4* 

tion  is  generated  by  a^a^  for  i,j  nearest  neighbours. 

We  look  at  the  form  the  kinetic  energy  takes  in  the 

continuum 


ft 


di p 


t 


9x 


9_£_ 

9  x 


dx 


K.E  . 


2m 


/ 


(2.2.13) 
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For  free  particles 


.  /  x  1  ^  -iKr 

ip  (r)  =  —  Z  a  e  — 

/V 


(2.2.14) 


therefore 


K.E  = 


2 ,  ,2 

fi  |  K[ 

2m 


Z  aKaK 


(2.2. 15) 


if  we  define  for  the  lattice  fluid  the  operators 

iKr  . 


aK  = 


t 

aK  = 


/N 

1 


Z  e 

j 


a  . 
3 


(2.2.16) 


/N  j 


-iKr .  . 

Z  e  a+ 

j 


One  possible  kinetic  energy  operator  with  the  correct 
continuum  limit  would  be 


K.E  = 


N-h2d 

2  . 

2ma  q  <  ij> 


(at  -  at )  (a  .  -  a  .  ) 

i  i  i  1 


(2.2.17) 


where  we  replaced  by  differences,  d  is  the  dimen- 

oX 

sionality,  q  the  number  of  nearest  neighbours 

2  iKr .  iKr .  , 

—  Z  (Z  (e  1  -  e  D )  aT)  x 

K 


K.E  = 


ft 


2 


2ma  q  <ij>  K 


-iKr  . 

x  (Z  (e  1  -  e 

K' 


-iKr  . 
—3 


)aK,) 


fi2d 

2 

2ma  q 


f 

^  amam 1 
KK'  K  K  <ij> 


(e 


iKr  . 

--i 


-  e 


iKr  . 
—3 


)  x 


-iKr .  “iKr . 
x  (e  1  -  e 


) 


(2.2.18) 
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In  the  limit  | K  j  <<  —  for  ij  nearest  neighbours 


iK . r .  iK . r  . 
;  -  -1  _  e  "  ^ 


iK .  a  e 


iK.r  . 

-  -l 


(2.2. 19) 


therefore 


K.E 


it 


,  iKr .  -iK'r. 

1  s^a^,  E  ( K . a )  e  .e  (K.a) 


0  2  “  ,  “K~K  .  . 

2ma  KK  <i;j> 


fi  +  9 

- j  E  a 'a  E  (K.a) 

2ma  K  <  ij> 


For  a  continuum  limit  the  nearest  neighbours  are  isotro¬ 
pically  distributed,  so  E  (K.a)^  °c  (K)  ^a^  ^ 

<  ij> 


•h  t  i  2  t 

K.E  ~  K  E  a 'a 

2m  1  1  K  K 

i\ 


(2.2.20) 


The  quantum  lattice  gas  Hamiltonian  given  by  this  heu¬ 
ristic  arguments  is 


dh 


(  ~~2  )  Z  (ai“ 
2ma  q  < i j  > 


f  t 


a  .  )  (a  .  -  a  . ) 
3  i  J 


t  t 

v  E  a  .  a  .  a  .  a  . 
°  <ij>  1  1  3  3 


For  spin  % 


(2.2. 21) 


ata .  =  (ox  +  iaY)(ox  -  i  oY)=  1/2  (1-cr?)  , 

li  i  ii  i  i 


therefore 


■ 


. 
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fi2d 


q2ma‘ 


E 

<  i  j  > 


/Xx,  y  y. 
(a. a.  +  a . a . ) 
i  D  iD 


v 

o 


v  z 
E  a -cr 

<ij> 


j 


r 


(2.2.22) 

which  is  our  anisotropic  Heisenberg  model.  In  absence 
of  interactions  this  reduces  to  pure  XY. 

By  the  assumption  that  a  strong  XY  interaction 
compared  to  the  Ising  like  part  yields  a  behaviour  like  the 
pure  XY  model,  we  can  hope  that  a  pure  XY  model  will 
describe  a  quantum  fluid  with  weak  interaction  be  it 
repulsive  or  attractive. 

This  condition  that  |JjJ>|J||  I  i s  here 

2 

d  fi 

—  »  1  ,  (2.2.23) 

v  q  0  2 

o  2ma 


where  a  is  the  lattice  spacing,  d  is  the  dimensionality, 
q  is  the  number  of  nearest  neighbours  and-vQ  the  poten¬ 
tial  energy  between  two  neighbouring  atoms. 

The  phase  diagram  for  the  superfluid  that  corres¬ 
ponds  to  this  hamiltonian  (2.2.22)  can  be  formed  in  the 
same  way  as  the  classical  gas  diagram  was  formed  in 
analogy  to  the  uniaxial  Ising  magnet 


/N 

1 

T 1 

s 

4/ 

U. 

T c  T* 

' 


■ 
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Our  hamiltonian  has  a  positive / ferromagnetic, 
J  as  it  came  from  the  kinetic  energy  and  is 


J 


q2ma 


2 


The  phase  diagrams  are 


\ 


The  pressure  in  the  fluid  still  corresponds  to  a  field 
in  the  z  direction  as  it  is  the  quantity  conjugate  to 

4*  z 

the  densityy which  is  E  aj_ai  ~  ^  °  •  Applying  a  field 
in  the  z  direction  would  not  destroy  the  order  until 
the  field  become  large  enough  to  rotate  the  magnetiza¬ 
tion  to  the  z  axis.  The  magnetic  fields  in  the  x  or  y 
direction  correspond  to  off  diagonal  fields,  \m  intro¬ 
duced  by  Bogolyubov  (1960) ,  Hohenberg  and  Martin  (1965) 
which  have  no  physical  realization. 


v^r 

o 


E  (v  •  a  .  +  v  .  a  .  ) 
.ii  ii 

l 


/ 


(2.2.24) 
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then 


9  v*  Nvo 


kT  n  1 

log  Z  =  — —  <aQ>  =  <  \jj> 


AT 


(2.2.25) 


o 


The  superfluid  order  parameter  can  be  defined  by 


|^(T,y)|  =  lim  |^(T,y,v*,v)|  .  (2.2.26) 

v ,  v*->0 

This  is  analogous  to  lim  <  M  > 

H  ->0 
x 
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2 . 3  Rigorous  one  dimensional  results 

Lieb/ Schultz  and  Mattis  (1961)  were  able  to  diagon¬ 
alize  the  Hamiltonian  of  the  anisotropic  XY  model  in  one 
dimension 

N 

21=  l  (1+Y)S*S*  +  ( 1  — y ) S^sY  i  (2.3.1) 

j  =  l  3  3+-L  3  3+-L 

by  using  the  transformation  to  fermion  operators  given 
by : 

at  =  exp  {ir  E  ctcv}  c^  ,  a|=ct  (2.3.2) 

3  K=1  R  R  3  11 

j=l  . 

a.  =  exp  {-iir  Z  V^}  e  ,  a^  .  (2.3.3) 

K —  1 

Z  Z 

They  calculated  the  correlation  functions  <  SqS~>  , 

X  X 

< Sq  S^>  and  proved  that  for. the  isotropic  case  .  y=0/ 
the  long  range  order  vanishes  even  at  T=0 

lim  <  SX  SX  >  ->  0  , 

o  2n 

n->°° 

while  it  does  not  for  y^O. 

Katsura  (1962)  independently  solved  the  same  model. 

z  z 

He  obtained  also  the  perpendicular  susceptibility  y  in 
finite  fields  and  the  ground  state.  Suzuki  (1966) 
obtained  the  staggered  initial  susceptibility  for  the 
isotropic  XY  model 


. 


. 


31 


X 


-  Nm' 
o  kT 


i/ 

7T 


dw 


o 


0  cosh  (K  cos  co) 


s  _  Nm' 
^o  ~  kT 


1  ^  tanh  (j(  cosw)  ^ 

tt  q  K  cos  o  daJ 


where 

K 


J 

kT 


(2.3.4) 


(2.3.5) 


(2.3.6) 


The  ground  state  functions  were  obtained  by  Katsura 


yo 

yo 


N/2 

(  n 

K>N/4 


A2K~1  A-2K+1) I 0> 


N/2 

(  n 

K>N/4 


A2K  A-2K)AnI0> 


for  ^  even  (2.3.7) 
for  |  odd  (2.3.8) 


where 

N 

A  =  —  Z  a  exp  [-i7T(^  -  ^)  ] 
K  /N  £=1  £  N  n  J 


1 

/N 


N 

E 

1=1 


+  r-  ' 
a£  exp  1_itt  (  — 


(2.3.9) 

(2.3.10) 


The  ground  state  energy  is 
E  0 

-Q-.-T  =  -  ±  .  (2.3.11) 

N  I  J  I  IT 

These  results  show  that  in  one  dimension  there  is  no 
long  range  order.  The  model  behaves  like  an  antiferro— 
magnet  with  respect  to  a  perpendicular  field  and  the 
ground  state  is  rather  like  that  of  the  anti ferromagnetic 
isotropic  Heisenberg  model. 


A  detailed  study  of  the  correlation  functions 
in  the  anisotropic  XY  model  was  done  by  Barouch  and 


32 


McCoy  (1969,  1970).  The  hamiltonian  treated  is: 

H=-J£[(l+Y>  S^S.^  +  (1-y)  SpY+1]-*y£S^  .  (2.3.12) 

They  studied  the  asymptotic  behaviour  of  the  correlations 

in  finite  fields  H  :  < SX  SX>  ,  < Sy  Sy>  ,  and  <SZ  SZ> 

z  o  R  oR  oR 

as  R  ->  °°  for  high  and  low  temperatures.  All  correlations 
decay  exponentially  with  the  exponent  depending  on  the 
field . 

For  T=0  the  behaviour  is  different  for  the  symme- 

2  2 

trie  cases  y=0,  or  'ftp  =J,  or  (-fty/J)  +  y  =1.  When 
none  of  those  conditions  obtain,  all  correlations 
approach  their  limit  exponentially  rapidly;  when  one  con¬ 
dition  holds  they  approach  the  limiting  values  as  some 
power  of  R.  The  effect  of  those  symmetries  can  be  seen 
at  finite  T  too  as  the  rate  of  the  vanishing  depends  on 
the  symmetry. 

The  parallel  susceptibility  has  not  been  calculated 
as  the  solution  of  the  XY  chain  with  a  field  in  the  z 
direction  used  all  simplifying  features  to  the  maximum,' 
the  nearest  neighbour  character  of  the  interaction,  the 
one  dimensionality  and  the  simple  quadratic  form  of  the 
Hamiltonian.  For  the  parallel  susceptibility  addition  of 
the  appropriate  Zeeman  term  gives  the  following  Hamiltonian 


, 

' 


I 


' 
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ft  =  -j  £  S*S*+1 


+ 


sysY 

bibi+l 


-  m  H  £  S 


x 


(2.3.13) 


Taking  x  ->  z ,  z  ->  y,  y  ->  x,  we  obtain 


ft =  -J  £  (a+a .  , ,  +  a . af , n )  -  J'  £  ata . af , , a .  , , +  J"  £  a+a . 

1  l+l  l  l+l  l  i  l+l  l+l  l  i 


(2.3.14) 


and  the  simple  quadratic  nature  of  ft,  is  destroyed. 


High  temperature  series  expansions  have  been  done 
by  Bonner  and  Fisher  (1964) .  They  calculated  the  energy , 
specific  heat.,  susceptibility  and  pair  correlations  for 
the  hamiltonian 


ji  =  -2J  lsisi+1  +  Y(S^S^+1  +  Spy+1)j  -  m  H  IS.  (2.3.15) 


for  Y=0  to  1,  therefore  not  reaching  the  region  of  XY 
symmetry . 


”  ■  :  :  '  ■  ’  '  £  '  ■:  i  -,5 
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2 . 4  Rigorous  thermodynamic  inequalities  and  scaling  laws 

Thermodynamic  arguments  were  used  by  Rushbrooke 
(1963,  1965)  and  Griffiths  (1965)  to  prove  rigorous 
inequalities  for  the  critical  exponents.  Most  of  those 
proofs  make  use  of  convexity  relations  (Ruelle  1963, 
Fisher  1964  and  Griffiths  1964)  of  the  free  energy 

F(aT^  +  £T2,  v)  ^  aF ( ,  v)  +  3F(T2,  v) 

F(T,  +  (3v2)  <c  aF(T,vx)  +  3F(T,  v2)  .  (2.4.1) 


The  Helmholtz  free  energy  A(T,v)  is  concave  in  the  tem¬ 
perature  and  convex  in  the  volume.  The  same  holds  with 
magnetization  instead  of  volume  for  a  hamiltonian  of 
the  form 


-  h  M  (2.4.2) 

or  provided  the  magnetization  commutes  with  the  hamil¬ 
tonian  (our  case  is  included  in  (2.4.2)),  we  will  write 
the  rest  in  magnetic  language. 

It  follows  from  the  convexity  properties  that 


(2.4.4) 


' 


' 


" 
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9S 
dTJ  H 


H 


therefore 


Since  C  ^  0 


C 

H 


T  ->  T  ,  H  =  0  . 
c 


(2.4.5) 


(2.4.6) 


Substituting  the  expected  behaviour  near  the  critical 
point  from  below  at  zero  field 


we  obtain 

f  (T)  (T  -  T)  ~a  >,  (T  -  t)2(^”1)+Y  g(T)  , 
o  c 

where  f  and  g  are  non-singular  at  T  .  Since  this  should 

close  to  T  it  follows  that 

n  / 


hold  for  all  T 


J 1 
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a'  +2(3+y'  >,2 


(2.4.7) 


This  is  usually  called  the  Rushbrooke  inequality. 

The  implicit  assumption  in  the  above  derivation 
was  the  existence  of  y',  that  is 

lim  Xm  T)  =  °°  for  T  Tc  . 

H+0 


The  same  assumption  and  use  of  convexity  of  A(T,  M) 
yield  (Griffiths  1965) 


a 


+  3(1  +  6)  >2 


(2.4.8) 


Many  other  inequalities  were  obtained  such  as 


y*  >>  3(6  -  1) 


T 


(6  +  1)  ^  (2  -  a)  (6-1)  • 


(2.4.9) 


(2.4.10) 


The  assumptions  made  to  obtain  them  are  plausible  for 
most  systems/such  as  a'4  ct  or 

( 9M^-)  <  0  for  H  ^  0 

9  H 

Buckingham  and  Gunton  (1969)  and  Fisher  (1969)  have 
proved  the  folio winq  inequalities 


(2  -  n )  v  >-  y 


(2.4.11) 


d 


6  -  1 
6  +  1 


>„  2 


n 


(2.4.12) 


■ 


, 


and 
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2^  2  -  n  (2.-4a3) 

where  d  is  the  dimensionality  of  the  system. 

The  additional  assumptions  that  went  into  those 
two  inequalities  1 2 ,  13  are  that  for  all  temperatures 
and  H  ^  0 


r2(r)  -  T*  >  Q 


r,  (r 


1 


r  3 }  " 


r. 


(ri) 


r. 


(r3" 


>>  0 


(2.4.14) 


where  T2  is  the  2  spin  correlation  function  and  therefore 
1^  is  proportional  to  the  spontaneous  magnetization, 

T1(TfE),  T2(T,H)  and  T4(T#H)  (2.4.15) 

are  monotonic  non-decreasing  functions  of  H, 


ri(T,H),  r2(T,H)  and  T4(T,H)  (2.4.16) 

are  monotonic  non— increasing  functions  of  T  for  any  field 
H  >>  0 .  These  assumptions  have  been  proved  for  the  Ising 
model  (Griffiths  (1967,  1970)). 

The  Josephson  inequalities  (1967)  are 


' 
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dv'  >  2  -  a' 


(2.4.17) 


dv  >  2  -  a 


His  assumptions  are  very  general  such  as  constant  volume 

or  constant  pressure  constraints. 

One  way  of  stating  the  basic  assumption  of  static 

scaling  theory  (Griffiths  (1968))  is  to  require  that  the 

singular  part  of  the  Helmholtz  potential  A  be  a  homogeneous 

1/8 

function  of  e  and  |m| 

A(As|xeM)  =  APA(s,  M)  .  (2.4.18) 

p  is  an  arbitrary  scaling  parameter  and  e  =(T-Tc)/Tc. 

This  means  that,  over  the  M,s  planey the  singular  part 
of  A  in  an  area  around  the  origin  will  be  the  same  as  for 
a  larger  scaled  area  when  itself  properly  scaled.  This 
was  first  proposed  by  Widom  (1965)  and  given  physical 
meaning  by  Kadanoff  (1966)  who  derived  this  homogeneity 
for  the  Ising  model  by  dividing  the  system  into  cells  of 
L  lattice  sites  per  side,  with  L  large  but  smaller  than  ~  * 

the  coherence  length,  a  is  the  lattice  spacing) ,  obtaining  the 
free  energy  and  claiming  that  it  should  not  depend  on  L. 

Taking  first  derivative  of  (2.4.18)  with  respect 

to  M  and  substituting 


■ 


39 


we  have 


A3  H  ( A  e  ,  A3M)  =  AP  H(e,  M) 


(2.4.19) 


Taking  a  second  derivative  with  respect  to  M  and  using 

i_  =  ,3H,  .  .  . 

a  3M.  T  f  obtain 

X>p 


A3  A3  M)  =  AP  M)  .  (2.4.20) 

On  the  critical  isotherm  H  ~  M  ,  setting  that  in  (2.1.20) 
equating  powers  of  A 

3  +  35  =  p 

Setting  M=0  and  T  ->  in  (2.1.18)  yields 

2  3  +  y'  =  p  .  (2.4.21) 

For  M=0,  T  T+  we  obtain  2  3  +Y=  p.  Eliminating  p  gives 

c* 

the  Widom  relation 


y  '  =  3(6  -  1 )  (2.4.22) 

and 

y  =  y'  .  (2.4.23) 

Taking  more  derivatives  yields  more  such  relations. 

All  the  rigorous  inequalities  we  quoted  in  the  first 
part  of  this  section  become  equalities  under  the  scaling 
hypothesis  (there  are  others  we  did  not  mention  that  do 
not) .  Since  A  is  a  function  of  two  variables  we  can  see 


' 
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that  only  two  independent  critical  exponents  are  allowed 
by  the  scaling  hypothesis. 

Taking  the  derivatives  of  the  free  energy  and  using 
the  same  homogeneity  argument  one  has  a  relation  for  the 
gap  exponents  defined  in  (2.1.7) 

F(3) (e)  =  x2p_33F(3) (As) 


but 


2  ,  . 
F  (e) 


-a'-y' 

(-£) 


SO 

=  2p  -  33  -  y'  =  3  +  Y ' 


but 

A ^  =  2  -  a '  -  3  =  3  +  Y' 

and 

A^  =  3  +  Y' 


Similarly  it  can  be  shown  that  all  gap  parameters  for  T< 

and  all  even  ones  for  T  >  T  are  equal. 

c 

These  results  of  scaling  hold  exactly  in  the  2- 
dimensional  Ising  model,  and  are  approximately  true  for 
3 -dimensional  Ising  and  Heisenberg  models.  Some  dis¬ 
crepancies  can  be  explained  by  allowing  larger  errors 
in  the  numerical  calculations. 
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The  Kadanoff  construction  applied  to  correlation 
functions  yields  a  scaling  behaviour  for  them  which 
already  for  the  3-dimensional  Ising  model  is  outside 
the  range  predicted  by  numerical  calculations. 

Since  these  arguments  are  extendable  to  dynamical 
scaling  we  will  deal  with  the  correlation  scaling  in 
the  dynamical  section  of  this  thesis,  and  mention  here 
only  that  the  equalities  obtained  for  the  static  indices 
are  (2.4.11)  and  (2.4.17)  with  the  equality  holding, 
that  is: 


dv '  =  2  -  a ' 

(2-n) v '  =  y ' 

d  =  Jbd. 

6  +  1  2  3+Y ' 


dv  =  2  -  a 
( 2-g ) v  =  y 

-  n  . 


(2.4.24) 


When  we  assume  scaling  all  those  equalities  will  be  assumed. 


' 

* 
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CHAPTER  III 

HIGH  TEMPERATURE  SERIES  EXPANSIONS 


3 . 1  General  method  of  graphical  expansions 


The  free  energy  can  be  expressed  as  a  formal 

series  in  powers  of  T_1.  The  method  has  been  expounded 

by  Domb  (1960,1965)  and  by  Fisher  (1965)  in  their  review 

articles .  We  have  Z  —  Tr  exp  (  (— It/kT)  and 

-  JL  =  lim  ~  (log  Z) ,  therefore 
kT  „ . _  N 


F 

kT 


lim  h  log  A  h  Tr  pn 

N->°°  n=0 


(3.1.1) 


where  K  =  J/kT  and  P  is  defined  by  K  =  - JP •  As  shown  by 
Rushbrooke  and  Wood  (1958)  the  logarithm  has  the  effect 
of  taking  only  the  terms  linear  in  N.  The  free  energy 
being  an  extensive  variable  this  was  to  be  expected. 

°°  Kn 

log  Z  =  2  log  N  +  log  (1  +  l  Vn  (3.1.2 


Tr  Pn 

where 


N 


For  the  spin  A  case  Tr  I  -  2  ,  I 


being  the  direct  product  of  N  spinor  unit  matrices. 

The  coefficients  un  are  polynomials  in  N  and  we  assume 
that  the  logarithm  is  also  expandable  in  a  power  series 


in  K 


oo  ^n 

log  (1  +  £  yn  ppr) 

n=l 


00 

E 

n=l 


x 

n  n . 


(3.1.3) 


i 
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which  is  a  definition  of  the  coefficients  A. 
result  can  be  checked  using 


n 


The 


n  ,  -i  xm+l  ,  > 
,  v  -1  m 

A  =  £  -  y 

n  ,  m  n 

m=l 


but  has  been  proved  generally.  Hence 


oo  n 

—  =  E  ^-r  Tr*  Pn 
kT  n=0  n* 


(3.1.4) 


.  m  (part  linear  m  N  of  Tr  A) 

where  Tr*A  =  - - - = - 

Tr  I 


Our  hamiltonian,  like  all  hamiltonians  used  in 
critical  phenomena,  is  a  sum  of  pairwise  interactions; 
therefore  each  term  in  Tr*Pn  can  be  represented  diagram- 
matically  as  a  linear  graph  of  n  lines  and  v  2n 

vertices . 

Tr*Pn  =  Z  w  (g^)  (g*?  ,  L)  .  (3.1.5) 

i 

g*?  are  all  possible  linear  graphs  of  n  lines  (i  is  a 
dummy  index  ordering  them)  .  (g_^,  L)  is  the  lattice 

constant  of  the  graph  g  determined  by  the  number  of 
ways  the  graph  gn  can  be  embedded  in  a  lattice  L. 
w(gn)  is  a  weight  assigned  to  the  graph, independently 
of  the  lattice , depending  on  the  interaction  and  the 


graph . 


' 
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3 • 2  The  graphical  expansion  method  used  for  the  XY  model 


The  properties  of  the  operators  a±  and  a^  simplify 
the  calculations.  The  method  used  by  Betts,  Elliott  and 
Lee  (1970)  is  as  follows:  The  general  term  wanted  for 
the  free  energy  is  Tr*  Pn.  Each  such  term  will  be  repre¬ 
sented  by  graphs  with  n  arrows  as  each  power  of  P  contri- 
butes  a  factor  a^a^  which  can  be  seen  as  an  arrow  from 
i  to  j .  Therefore  a  priori  all  linear  graphs  of  n  arrows 
and  v  vertices  v^2n,  can  be  expected  to  appear  in  Tr*  Pn. 
A  non  zero  contribution  can  only  be  obtained  when: 

(a)  an  equal  number  of  arrows  enter  and  leave  each  site; 

(b)  arrow  heads  and  tails  alternate  at  each  site. 

This  follows  because 


a. a.  =  a tat  =0  (3.2.1) 

i  i  i  i 

and 

Tr  a .  =  Tr  a !  =  0  .  (3.2.2) 

l  l 

Were  there  an  odd  number  of  operators  at  a  site  we  could 
reduce  them,  two  at  a  time  using  {a^,  a^}  =  1,  to  a 
single  traceless  operator. 

The  partition  function  is  the  sum  of  all  ordered 
directed  graphs  satisfying  (a) .  Several  of  those  will 
correspond  to  one  directed  shadow  graph.  Several  shadow 
graphs  will  correspond  to  one  bare  graph  in  which  bonds 
replace  the  directed  arrows. 


' 
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Example 


the  shadow  graph 


Out  of  51=120  possible  orderings  of  the  5  arrows  only  20 
will  give  a  non  zero  contribution  and  they  all  correspond 
to  the  same  shadow  graph. 

This  shadow  graph  is  equivalent  to  5  other  shadow 
graphs,  the  following  two 


2 


2 


1 


and  those  obtained  from  these  three  by  reversing  all  arrows. 
All  six  are  represented  by  the  same  bare  graph 


2 


1 


3 


The  trace  of  each  such  graph  is  2N  S,  where  s  is 


the  number  of  vertices  of  the  bare  graph,  since  for 

t 

spin  h  the  trace  of  the  unit  matrix  is  2  whrle  Tr  a^^ 


( 
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4* 

Tr  a^a^  -  1*  This  is  a  great  simplification  compared 
to  the  Heisenberg  model  where  traces  of  different  spin 
operators  have  to  be  calculated. 

The  horizontal  weight  of  each  shadow  graph  g' 
is  defined  by 

h(g')  =  e(g')  S(g)/S(g')  (3.2.3) 

where  S(g)  is  the  symmetry  number  of  the  bare  graph  g, 
S(g')  is  the  symmetry  number  of  the  shadow  directed 
graph  g'  and  e(g')  equals  one  if  g'  is  equivalent  to 
the  graph  obtained  from  it  by  reversing  all  arrows. 

Going  back  to  the  above  example,  the  symmetry  number 
of  the  triangle  is  6,  and  the  symmetry  number  of  our 
graph  is  1.  e(g')  =  1  as  in  this  case  interchange  of 

two  labels  is  the  same  as  reversal  of  arrows.  Equivalence 
of  shadow  graphs  is  defined  as  follows:  two  graphs  are 
equivalent  if  their  arrows  and  vertices  are  in  one  to  one 
correspondence.  Such  a  correspondence  is  an  embedding. 

The  number  of  embeddings  of  one  graph  in  itself  is  its 
symmetry  number . 

The  vertical  weight  v (g1 )  is  defined  as  the  number 
of  allowed  ordered  directed  graphs  corresponding  to  the 
shadow  graph  g' .  In  our  example  that  number  was  20. 
Allowed  here  meant  giving  a  non  zero  contribution,  there¬ 
fore  satisfying  criterion  (b) . 
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All  these  weights  are  calculated  by  computer;  for 
that  purpose  graphs  are  represented,  as  usual,  by  ma¬ 
trices  or  vectors.  For  the  directed  shadow  graph  the 
matrix  element  M(i,j)  is  the  number  of  arrows  from  i  to 
j .  The  bare  graph  matrix  has  a  unit  element  in  a  location 
i,j  if  the  vertices  i  and  j  are  connected,  and  zero  other¬ 
wise.  For  our  example 


'0 

1 

1' 

'0 

1 

1' 

M(g')  = 

0 

0 

1 

,  M(g)  = 

1 

0 

1 

2 

0 

0 

> 

1 

1 

0 

In  order  to  exclude  equivalent  graphs,  to  find  among 
a  group  of  lattice  constants  the  one  pertaining  to  the 
bare  graph  and  to  calculate  the  horizontal  weight  the 
incidence  matrices  defined  above  were  used. 

For  the  vertical  weight  another  kind  of  matrix 
representation  was  used.  The  matrix  is  rectangular. 

Rows  denote  order  of  arrows,  and  the  columns  specify 
vertices.  The  arrow  in  the  k ' th  position  being  i-j  the 
matrix  elements  in  the  k 1 th  row  will  be  N(k,i)=l,  N(k,j)=-1. 
The  vertical  weight  of  the  graph  corresponds  therefore  to 
the  number  of  permutations  of  rows  of  the  matrix  N  that 
alternate  1  and  -1  in  each  column.  Full  use  of  efficient 
methods  of  counting  permutations  was  made  in  the 
computer  program  written  by  Dr.  C.J.  Elliott  for  the 
free  energy  calculation.  This  calculation  being  rather 


■ 
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time  consuming  for  higher  order  graphs,  several  theorems 
were  proved  (Betts,  Elliott  and  Lee  (1970))  which  reduced 
noticeably  the  computer  time  required  and  were  based  on 
full  use  of  the  fact  that  there  are  no  restrictions  at 
vertices  of  order  2 .  We  will  mention  those  used  and 
their  generalizations  for  our  calculations  in  the  next 
section . 

To  summarize/ the  final  contribution  to  the  term 
Tr*  Pn  is 


Z 


(g  .  , L) v  (  g '  •  ) h ( g '  .)  2 


-s 


n,  i 


' n ,  i 


(3.2.4) 


gn,i 

where  the  sum  is  over  all  inequivalent  directed  shadow 
graphs  of  n  arrows/  g^  ^ .  Note  that  the  lattice  constant 
is  that  of  the  bare  graph  g  .  the  horizontal  weight 

1*1  f  1 

giving  the  occurrence  of  the. directed  graph  g^  ^  for 
each  bare  graph  ,  the  vertical  weight  being  the  number  of 

—  g 

allowed  orders  of  the  arrows  and  2  being  the  actual 
trace . 


■ 
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3 • 3  The  series  expansion  method  for  the  initial  parallel 

susceptibility 

A.  The  XY  model  is  invariant  under  rotation  in  the  xy 
plane,  therefore  we  can  choose  any  direction  in  that  plane 
to  calculate  the  initial  parallel  susceptibility.  We  add 
a  Zeeman  term  and  the  resulting  hamiltonian  is: 


$i  =  -  J  I  (a+a  .  +  a.a+)  -  H  M 

<ij>  1  ^  ID 


(3.3.1) 


The  initial  parallel  susceptibility  is  defined  by 


o  _  xx  _  jL_ 
Al  ^o  N3 


3H 


2  log  Z 


x 


(3.3.2) 


H  =0 
x 


r  -i  12 

For  models  where  ~  0  this  reduces  to  —  <M^>  r  which 

is  the  square  of  the  fluctuation  in  the  long  range  order. 
The  relation  between  the  two  is  discussed  in  Subsection 
3.3B.  Using  the  notation  of  the  last  section,  Jt  =  -  J  P , 
we  have : 


X, 


o 


1  3' 


NB  3h3 

X 


log  Tr  e 


3 (JP  +  H  M  ) 
xx 


H  =0 
x 


(3.3.3) 


As  before  we  use  the  relation 


1  .  m  3A 
—  log  Tr  e 
N 


,n 


=  Z  Tr*  An 
n ; 


(3.3.4) 
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and  obtain: 


X 


o 


1  3 

3 


9H 


n 

T  I  t  Tr* (JP  +  H  M  ) n 
^  n:  xx 

x 


H  =0 
x 


(3.3.5) 


We  take  the  second  derivative  and  set  H  =0.  This 

x 

means  keeping  only  the  term  proportional  to  H2 .  The 

2 

coefficient  of  H  in  Tr*  (JP  +  H  M  ) n  can  be  written 

x  xx 

in  the  following  form: 

n-2  ,  _?_v 

h  Tr*  I  n  ( JP )  M  (JP)n  M  .  (3.3.6) 

k=0  x  x 


We  obtain  (3.3.6)  as  follows:  The  term  with  two  factors 
B  in  Tr  (A  +  B)n  is 

n-s-2  n-2  „  0  n 

Tr  Z  Z  A£  BAn  2  1  S  BAS 

£=0  s=0 

We  use  the  cyclic  property  of  the  trace  to  bring  the  factor 
A  to  the  front.  The  last  expression  reduces  to 

n-2  .  9 _ • 

Tr  Z  (j+1)  AD  BA  2  3  B  (3.3.7) 

j  =  0 

as  j  =  £+s  can  be  obtained  in  j+1  ways  by  adding  two 
integers  £  and  s.  Again  we  invoke  the  cyclic  property 

*  2  *  *  i  2  *  i 

of  the  trace  and  note  that  Tr  A-^  B  An  ^  B  =  Tr  A2  B  An  ^  B 
for  j'=  n-2-j .  We  add  those  two  terms  and  divide  by  two. 


. 
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%[<j  +  l>  +  (n-2-j+l)]  Tr  A-^  BAn  2  ^  B 

=  %  n  Tr  Aj  BAn_2"3  B  .  (3.3.8) 

For  even  n  the  sum  (3.3.7)  has  an  odd  number  of  terms. 
Still  the  expression  (3.3.8)  holds  because  the  middle 
term  occurs  for  j  =  and  is  j  Tr  A^  BAn~2"^  B. 

Now  that  we  have  justified  the  use  of  (3.3.6)  we 
can  substitute  it  in  (3.3.5)  and  obtain 


X°  =  £  E  ggp  ni2  Tr*  Pk  M  p-2-k  M 
J  n  <n  •  k-0  x  x 


(3.3.9) 


m  f 

We  have  M  =  ^  E  (a!  +  a.)  ,  hence 

X  Z.  11 


X°  =  j  £  §t-  f  .E.  Tr*(Pka+Pn  1  ka.+  Pka . Pn  1  kaf) 
J  n  n*  k=0  1  3  i  3 


(3.3.10) 

*4*  XI.  “  1 _ k 

Each  term  in  Tr*  (P  al  a^)  can  be  represented  by  a 

graph  with  n-1  solid  arrows  and  a  dotted  arrow.  The  dotted 

f 

arrow  represents  the  term  a. a.  which  is  not  restricted  to 

1  3 

nearest  neighbours.  For  i=j  we  have  a  graph  of  n-1  solid 
arrows  and  a  dotted  "bubble"  which  can  be  on  a  vertex  of  the 


graph  or  off  the  graph. 

The  calculation  is  done  separately  for  the  three 


cases : 
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(a)  i^j  Each  ij  term  occurs  twice  therefore  the  extra 

3.  r_ — — ,3  factor  2  in  (3.3.12). 


4v 


I  L - <(. - 4 

(b)  i=j  and 


f  f  t  t 
ala2a2a3a4a3a4al 


i  is  on  the  graph 


t  t  t  f 

a1a2a3a2a3a3a1a3 


i  is  not  on  the  graph 


t  t  t  f 

al  2a2  3a3ala4a4 


Case  (c)  is  the  simplest.  Since  i  is  not  on  the  graph 


9n-l 


+ 

,  a^  and  a^  commute  with  the  factors  P  and 


m  *  /T,k  t„n-l-k  ,  _k  n-l-k  t.  m  *  „n-l 
Tr*  (Pa.P  a.+Pa.P  a.)  =  Tr*  P 

l  li  j 


But  here  Tr*  Pn  ^  does  not  give  the  same  contribution 
as  5fn_2_  gave  in  the  partition  function.  The  lattice 
constant  is  now  the  separated  lattice  constant  (g^_^.;L), 
which  is  ~s(gn_-)  ;L)  .  s  is  the  number  of  vertices  of 


g  i  • 

n-l 

The  other  factors,  horizontal  and  vertical  weights 
and  trace  are  the  same  as  in  the  partition  function. 

The  resulting  series  is: 


X 


o  (c) 

II 


K11'1 

(n-l)  ! 


E  ■ 

9n-l 


h(gn-l)(gn-l;L)s  V(gn-1> 


(3.3.11) 


The  weights  are  given  in  Table  A. 9  in  the  Appendix. 
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For  case  (a)  i^j ,  the  dotted  arrow  is  placed  with 
the  head  at  the  last  level  and  we  sum  over  all  allowed 
levels  for  its  tail.  Since  vertices  of  order  2  add  no 
restrictions  we  introduce,  for  the  vertical  weight 
calculation  only,  a  spurious  vertex  in  the  dotted  arrow. 


+  f  f  f  t 
ala2a2a3a4a5a3a4a5al 


We  allow  the  new  arrow  formed  by  the  former  tail  and  the 
spurious  vertex  to  permute  freely  with  the  other  solid 
arrows.  The  arrow  formed  by  the  spurious  vertex  and  the 
tail  is  kept  fixed  at  the  last  level. 

The  lattice  constant  is  calculated  after  deleting 
the  dotted  arrow  since  it  is  not  restricted  to  nearest 
neighbours.  The  horizontal  weight  is  as  defined  in 
(3.2.3).  In  order  to  keep  the  dotted  arrow  different 
from  the  others  we  put  a  minus  sign  in  front  of  the 
matrix  element  M(i,j).  N(i,j)  and  N(j,i)  become  zero 
if  no  other  arrow  connects  i  and  j . 

The  graphs  were  obtained  by  Betts,  Elliott  and  Lee 
(1970)  for  the  fluctuation  series.  We  calculated  the 
new  vertical  weights  only 


X 


o  (a) 


2 

J 


n 


K_ 

n 


n 


h(gn-l)v(gn-l)  (gn-l'L) 


(3.3.12) 


n-1 


' 


.* 
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where  9n_^  is  the  graph  with  n-1  solid  arrows  and  a 
dotted  arrow.  gn-1  is  the  bare  graph  with  the  dotted 
arrow  removed.  The  weights  are  given  in  Table  A. 5. 

In  case  (b)  we  have  i =  j ,  i  eg.  The  same  pro¬ 
cedure  as  for  case  (a)  is  followed  for  the  vertical 
weight  calculation.  Lattice  constants  are  those  of  the 
bare  graph  without  the  bubble.  The  graphs  were  generated 
by  taking  the  partition  function  graph  adding  a  "bubble" 
at  each  vertex  and  eliminating  equivalent  graphs.  We 
kept  the  spurious  vertex,  s+1,  in  the  generation  process 
and  the  horizontal  weight  calculation.  In  order  to  keep 
the  bubble  different  we  introduced  a  minus  sign  in 
M ( i , s+1 ) ,  M(s+l,i)  and  set  N (i , s+1) =N (s+1 , i) =0 ,  that 
is  the  bubble  is  deleted  in  the  bare  graph.  The  result 
is  reflected  in  the  following  example 


Let  g^_^  be  a  partition  function  graph  9n_-^  with  a 
bubble  on,  then: 


X 


o  (b) 


1 

J  L  n! 


h  (gb  ,  )  v  (gb  ,  )  (g 
^n-1  ^n-1  ^n- 


•1 


,L) 


gb  , 
^n-1 


(3.3.13) 


The  weights  can  be  calculated  from  Table  A. 7. 
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The  calculation  of  fluctuation  type  vertical  weights 
(one  arrow  fixed)  is  facilitated  by  the  lemmas  and  theorems 
of  Betts,  Elliott  and  Lee  (1970)  . 

Lemma  1 

Let  g  be  a  graph  with  £  arrows  and  g*  is  the  graph 
obtained  from  g  by  turning  one  arrow  into  a  dotted  arrow. 
Let  there  be  £^  arrows  from  vertex  a  to  vertex  b  one  of 
which  is  the  dotted  arrow.  Then 

_  Vt 

vg*  i 


Proof 

If  £^=  1  fixing  one  arrow  at  a  level  reduces  the 
number  of  possible  vertical  orderings  by  y  .  The  reason 
is  in  the  cyclical  property  of  the  trace,  all  cyclic 
permutations  of  an  allowed  permutation  are  allowed. 

If  £^  >  1  there  are  £^  arrows  which  were  indistinguishable 
in  g  and  by  fixing  one  we  made  them  distinguishable. 

Lemma  2 

Let  g  be  a  graph  with  £  arrows  and  s  vertices  and 
let  one  arrow  be  of  order  2  in  both  vertices.  Then,  one 
vertex  of  order  2  can  be  eliminated.  The  resulting 
graph  of  £-1  arrows  and  s-1  vertices  has  a  vertical 
weight  v^ . 

=  v^ .  £ 


v 


, 


56 


The  procedure  can  be  repeated  and  all  arrows  of  order  2 
in  both  vertices  eliminated.  The  proof  is  obvious. 

Theorem  I 

Let  g  be  a  disconnected  shadow  graph  of  £  arrows 
consisting  of  two  subgraphs  g^  and  g^  having  no  vertices 
in  common.  If  g^  and  g^f  consisting  of  £^  and  £2  arrows 
respectively,  have  vertical  weights  v^  and  ,  then  the 
vertical  weight,  v,  of  g  is  given  by 

v  =  (£:/£1!  £2'.  )  v x  v2  .  (3.3.14) 


Proof 

As  there  is  no  interference  between  the  arrangements 
in  the  two  subgraphs  the  vertical  weight  of  the  graph  g 
is  given  by  the  product  of  two  factors.  The  first  is  the 
number  of  ways  of  assigning  £  vertical  levels  such  that 
£-^  are  occupied  by  the  arrows  of  g^  and  the  remaining 
£2  =  £  -  £-^  are  occupied  by  the  arrows  of  The  other 

factor  is  then  the  product  of  the  vertical  weights  of  the 
component  subgraphs. 

Theorem  II 

Let  a  connected  shadow  graph  g  of  £  arrows  be 
composed  of  two  subgraphs  g^  and  g2  of  £^  and  £2  arrows 
respectively  such  that  g^  and  g2  have  m  common  vertices 
of  order  2  in  g  (hence  of  order  1  in  each  of  g1  and  g  ) 


■ 
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and  no  other  common  vertices.  Then  the  vertical  weight, 
v,  of  g  is  given  by  the  same  expression  as  in  theorem  I. 

Proof 

As  there  is  no  restriction  on  head  and  tail  vertical 
sequence  at  vertices  of  order  2,  the  vertical  weight  of  g 
must  be  the  same  as  that  of  a  separated  graph  g'  consis¬ 
ting  of  subgraphs  g^  and  g^.  But  the  vertical  weight  of 
g  is  given  by  (3.3.14)  and  so  the  theorem  is  established. 

Theorem  III 

Consider  a  shadow  graph  g  of  order  £  composed  of 
two  subgraphs  g^  and  g^  of  £^  -  1  and  £^  -1  arrows 
respectively  and  a  third  subgraph  g3  consisting  of  a 
single  arrow.  The  subgraphs  g^  and  g^  have  in  common 
m  vertices  of  order  2  in  g  and  no  other  common  vertices. 

The  subgraph  g^  (the  single  arrow)  connects  vertices  of 
arbitrary  (even)  order  orders  r^  and  r2  in  g^  and  g2 
respectively.  Then  the  vertical  weight  of  g  is  given 
by  (3.3.14)  again.  However,  now  v^  and  v2  are  the  vertical 
weights  of  the  graphs  gj  and  g^  where  gj  is  a  subgraph  of 
g  consisting  of  g^  and  g^ ,  and  g!,  is  a  subgraph  consisting 
of  g^  and  g^ .  It  is  important  also  to  note  that  £^  +  = 

£  +  1  . 

Proof 

Observe  first  that  the  vertical  weight  of  g  is  the 
same  as  that  of  the  graph  g'  obtained  from  g  by  separating 


' 
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all  the  second  order  vertices.  Next  observe  that  in  counting 
the  number  of  allowed  arrow  configurations  for  any  shadow 
graph  of  order  Z  we  may  assume  that  one  arrow  is  at  a  fixed 
level  and  compute  the  number,  n,  of  allowed  configurations 
of  the  remaining  arrows.  Then  the  vertical  weight,  v=£n. 

In  g'  let  the  arrow  forming  g^  be  fixed.  The  remain¬ 
ing  Z  -  1  arrows  belong  either  to  gx  or  g2«  The  £  -1 
arrows  may  be  arranged  among  the  Z ^  -1  levels  assigned 
to  g-^  in  n^  ways  independently  of  the  arrangement  in  n2 
ways  of  the  Z^  -  1  arrows  of  g2  among  the  -  i  levels 
assigned  to  g2 .  The  Z  -  1  levels  may  be  assigned  in 
(Z  -1  )  1/ -1  )  !  (&2  -I  )  ways  giving 


n  =  nln2  '  (^2_1) 


(3.3.15) 


Noting  that  n^  =  v^/£^  and  n2=  v2/£2  the  result  (3.3.14) 
follows . 

B.  We  discuss  here  the  relation  of  the  fluctuation  to 
the  intial  susceptibility. 

As  seen  in  Sec .  5.4 


(3.3.16) 


0 


but 


+  [#,M]]  +  ...  (3.3.17) 


=  M  +  X  fy,  m] 


. 

jn 
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therefore 


X 


T 


3  Y  + 


<  M  [1C,M]>  + 


(3.3.18) 


Each  successive  term  in  (3.3.18),  when  expanded  in  a  high 
temperature  series,  starts  its  contribution  at  a  higher 
power  of  3  and  its  contribution  is  less  divergent.  The 
reason  is  that  in  a  term  <  M[J{,  [JC,  [.  .  .  ]]...]  >  each 

commutator  has  the  effect  of  removing  two  spins. 

Falk  and  Bruch  (1969)  proved  that  for  a  general 
thermodynamical  system 


-k 


8oo 


1  -  e 


X 


T 


8  a) 


3Y 


«  1 


where 


(3.3.19) 


co 


<[m,  [#,m]]  > 

2  Y 


(3.3.20) 


The  equality  of  the  indices  for  Y  and  XT  follows  from 
the  assumption  that  the  numerator  in  oo  is  finite  as  T->  T  . 
Though  <  [m  ,  [jt  ,m]  J  >  has  been  shown  to  be  finite  for  the 
isotropic  Heisenberg  model  by  Mermin  and  Wagner  (1966), 
their  proof  cannot  be  generalized  to  the  parallel  suscep¬ 
tibility  of  the  XY  model. 

Our  results  indicate  that  <  [M ,  [# ,M]  _]  >  ,  which  is 
the  first  moment  of  the  frequency  dependent  susceptibility 


' 


' 
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is  not  finite.  As  T  -*  the  first  moment  diverges 

with  exponent  a~  ~  0.09.  Nevertheless  oo  — — *■  0 

2  T-*T 

c 

because  the  fluctuation,  Y,  diverges  much  more  strongly. 

The  proof  of  Falk  and  Bruch  follows  from  the 
following  expressions  for  the  susceptibility 

3/2  y  AO  _y.t l> 

XT=  2  /  dX  <e  J  M  e  M  >  (3.3.21) 

0 

and 

Xm  =  /  do)  S  (w)  -i— - - (3.3.22) 

T  -co  °  Jj  to 


where 

00 

SO)  =  JiN  /dt  e"1Ut  <M'M(t)> 

^  —00 

and 

oo 

/dm  S  (oo)  =  Y 

—  oo 

The  integrand  in  (3.3.21)  has  a  maximum  at  A=0,  from  which 
it  follows  that 

XT  <  Y 

The  second  inequality  follows  using  the  convexity  of 


x 
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3.4  Fourth  order  fluctuation 


The  gap  parameter  was  defined  in  (2.1.7).  To 
calculate  the  gap  parameter  one  has  to  calculate  the 
critical  index,  of 


F2  = 


1 


4  4 
Nm  3 


9  H 

x 


log  Z 


H  -0 
x 


(3.4.1) 


The  gap  parameter  A  is  then 


A  =  h  (y 2  ~  y) 


(3.4.2) 


We  note  that  upon  differentiation  (3.4.1)  becomes 


^4  3 (JP+M  H  ) 

— — r  Tr  e  x  x 


4  4 

Nm  3  F. 


9  H 

x 


H  =0 
x 


-  3 


Tr  e 


KP 


9H' 


Tr  e 


(JP+M  H  ) 
x  x 


x 


H  -0 
x 


Tr  e 


KP 


(3.4.3) 


To  obtain  (3.4.3)  we  made  use  of  the  fact  that  above  Tc 
and  at  H=0  odd  derivatives  of  Z  vanish. 

When  the  magnetization  and  Hamiltonian  commute, 

(3.4.3)  reduces  to 


34  {<m4>  -  2  <m2>2} 

x  x 


(3.4.4) 
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The 


second  term  in  (3.4.3)  is  proportional  to  (y^)^. 


In  Section  3 . 3B  we  saw  that  behaves  like  <M 


as  T  T  .  The  first  term  in  (3.4.3)  is 


d4  4  „n-4 

3  m  K 


,  ~T\~r  Tr  I  PDM  PkM  PVM  P  M 

(n—  4)  .  .  .  x  x  x  x 

n  ij  £k 


(3.4.5) 


The  difference  between  (3.4.5)  and  <M4>  is  in  terms  of 

x 

commutators  of  M  and  .  By  physical  reasoning  this 
difference  is  expected  not  to  diverge,  or  diverge  less 
than  the  same  quantity  with  the  commutator  brackets 
removed . 

As  mentioned  in  Section  3 . 3B  commutators  have 
the  effect  of  removing  two  spins  from  each  trace,  thereby- 
reducing  the  divergence. 

A  more  rigorous  treatment  means  going  to  the 
fourth  order  in  the  field  in  the  expression  for  the 
density  matrix 

4  t  M  t2  h  bp 

(jP  /  f  f  fe  |_Mx(t1)  ,  [Mx(t2)  ,  [Mx(t3)  ,  [Mx(t4)  ,  p]  ]  J  ] 

—  OO  —CO  — OO  —00 

-i  tJ[ 

x  e  ^  dt-^dt^dt^dt^  ,  (3.4.6) 


which  is  rather  complicated. 
We  treat  therefore 


1L.  i  {  1  <  k2>2 

4  N  1  2  x 
m 


^<M4>  } 

2  x 


(3.4.7) 


t 


The  definition  (3.4.7)  has  the  advantage  of  yielding 

a  positive  series  starting  with  unity. 

The  term  <M4>  is  a  sum  over  four  indices 
x 


<M4>  =  ~  Z  <CJ* 
x  16  i,  j  ,k,Jl 


a 


x 


a 


x 

k 


<t> 


We  consider  all  possible  cases  for  the  summand  and 
list  them  below: 


(i) 

(ii) 

(i-ii) 

(iv) 

(v) 


x4 

<o  > 
a 


^  x2  x2 
<  a  a,  >  , 

a  b 


x3  x. 

<a  a,  >  , 

a  b 


x2  x  x. 

<a  a,  a  >  , 

a  be 


x  x  x  x^ 

<  a  a,  a  a  ,>  , 

abed 


a^b 

a^b 

a^b^c 

a^b^c^d 


Summing  (i)  ~  (v)  separately,  we  obtain 


<M  *>  =  (N  +  3N  (N-l )  +  8N  Z  <  a*  a*> 


,4 .  m 


1  6 


j^O 


v  v  x  x  x  x  ■! 

+  12N (N-2 )  Z  <a  a.>  +  Z  <aiajak  £  } 

j^O  J  i^j^k^£ 


+ 


(3.4  . 


' 
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..2^  m  N  , ,  ,  „  „  xx, 
<M  >  =  ~r-  (1  +  2  £  <  a  o  ■>) 
x  4  j^O  °  2 


4  2 

3<M2>2  =  (3  +  12  E  <aXaX>  + 

x  16  j^O  °  ^ 


(3.4.9) 


.  ,  n  v  „  X  X^  „  X  X  . 

+  12  E  <  a  a.><a  o,  >) 
j^o  °  3  °  k 

k^O 


We  see  that  the  parts  of  the  second  and  fourth  terms  in 
(3.4.7)  which  were  proportional  to  N  are  cancelled  by 
the  first  two  in  (3.4.9).  The  third  term  of  (3.4.9) 

X  X  X  X 

subtracted  from  E  <a.a.a,  a„>  leaves  the  part  which 
is  linear  in  N. 

When  [M,  ,?f  ]  =0  there  is  no  need  to  check  on 
linearity  in  N,  as  is  a  derivative  of  the  free  energy 


which  we  know  to  be  linear  in  N. 

We  can  define  a  generalized  exponential  e 


A+B 


rn 


which  factors  B  will  always  be  taken  to  the  front 


A+B 


A  B 
e  .  e 


Using  this  definition  of  order  we  see  that 

.  4 


Y2  = 


16 


4  4 
Nm  3 


3  H 


n  m  -3  (X-M  H  ) 

—  log  Tr  e  x  x 

4  ^g  g 


x 


H  =0 
x 


Kubo  (1962)  discussed  these  generalized  exponentials 
and  the  calculation  of  generalized  cumulants  and  moments 
with  them.  The  relation  between  cumulants  and  moments 


■ 
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remains  unchanged  as  long  as  we  apply  the  generalized 
order  consistently -  Therefore 


,  m  -3  tft-M.  H 

log  Tr  e  x  x' 

g  g 


is  equal  to  the  cumulant  expansion  series  which  is  linear 
in  N. 

We  see  from  (3.4.8)  and  (3.4.9)  that 


= 


16m 

N 


{1  +  8 


<  o 


j^O 


x 

o 


aX> 

3 


-  Q  } 


(3.4.11) 


where 


Q  =  -h 


K 


n 


n  I 


m  .  n  x  x  x  x 
Tr*  P  a  .  a  .  a,  a  0 
l  3  k  £ 


The  second  term  in  (3.4.11)  is  four  times  the  flucutation 
series  apart  from  the  initial  unity.  We  calculate  the 
third  term.  We  generate  the  graphs  by  opening  two  bonds 
in  the  partition  function  graphs,  under  the  restriction 
that  the  four  vertices  be  distinct. 

I"1 

We  note  that  each  graph  P  aaaj-)aca(j  appears  24 

times  in  the  sum  as  each  index  i,  j,  k,  £  can  take  any 

value  a,  b,  c,  or  d.  On  the  other  hand  we  are  counting 

each  term  twice  in  the  graphical  expansion  because  the 

"f*  *t* 

dotted  arrow  associated  with  aaap>aca(3  can  be  either 


. 


' 


for  example 


a+b,  c-*-d  "  or  a->d ,  c-+b, 


Therefore 


Kn 

Q  =  -e  z  jri, 

n  gn+2 


h(gn+2)V(gn+2>  L> 


n 


(3.4.12) 


The  vertical  weight  is  calculated  with  the  dotted  arrows 
fixed  at  the  last  two  levels  (since  the  four  vertices 
are  distinct  their  order  is  immaterial) . 

For  elimination  of  equivalent  graphs  and  calcu 
lation  of  horizontal  weights  the  matrix  elements  M(a,b) , 
M(c,d)  are  preceded  by  a  minus  sign.  N(a,b) ,  N(b,a) , 
N(c,d)  and  N(d,c)  are  zero  if  the  vertices  a-b  and  c-d 
respectively  are  not  connected  by 


other  arrows . 


. 
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CHAPTER  IV 

THE  SERIES  FOR  THE  STATIC  QUANTITIES  AND 
THEIR  ANALYSIS 

4 . 1  Series  analysis 

The  functions  we  deal  with  are  assumed  to  have  near  Tc 
an  asymptotic  behaviour  of  the  form: 


(g+l) 


(4.1.1) 


f  (x)  A  (xc  "  x) 


We  have  a  finite  number  of  terms  in  a  series  expansion  of 
f(x)  in  powers  of  x.  Two  basic  methods  have  been  deve¬ 
loped  to  obtain  estimates  of  xc  and  g  from  this  truncated 
expansion 


f  (x)  =  Z  a  xn  +  R. 
n-0  n 


(4.1.2) 


N 


(a )  The  ratio  method 

This  method  was  first  used  by  Domb  and  Sykes  (1957)  . 
Had  (4.1.1)  been  an  exact  equality,  using  the  binomial 
expansion  we  would  see  that  the  ratios  r^  defined  by 


(4.1.3) 


would  be  equal  to 


1 


r 


n 


x 


c 


(4.1.4) 
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If  the  series  (4.1.2)  converges,  its  radius  of  conver¬ 
gence  in  the  complex  x  plane  is  the  distance  from  the 
origin  to  the  first  singularity.  The  limit  of  |  r  ^  |  < 
if  it  exists,  is  the  radius  of  convergence.  For  series 
with  all  positive  terms  the  first  singularity  must  lie 
on  the  positive  real  axis  and  can  be  expected  to  be  the 
physical  singularity  we  are  seeking;  therefore  we 

expect  the  ratios  to  approach  x  ^  as  n  ->  00 .  Since  we 

o 

expect  (4.1.1)  to  hold  near. x  the  relation  (4.1.4)  will 

obtain  in  the  limit  n  -*  00 .  A  plot  of  rn  against  1/n  will 

therefore  approach  a  straight  line  with  intercept  1/x 

at  i=0  and  a  limiting  slope  of  g. 
n 

Various  modifications  and  refinements  of  the  ratio 

method  have  been  introduced.  Neville  tables  have  been 

used  by  Baker,  Eve,  Gilbert  and  Rushbrooke  (1968)  which 

are  successive  linear  extrapolations.  The  quantities  ln, 

q  etc.  are  constructed  1  =  n  r  -  (n-1)  r  ,  , 

n  n  11  ± 

qn  =  %  [n  ln  -  (n-2)  q.J,  cn  =  i  [n  qn  -  (n-3)  g^J  , 
etc.,  and  as  long  as  they  are  monotonic  they  provide 
improved  estimates  of  x"1 .  Once  the  critical  point  is 
evaluated,  better  estimates  for  the  exponent  can  be 
obtained  by  the  relation  (Domb  and  Sykes  (1957)) 

n  txc(an/an-l)  _  d  +  g  as  n  -  «° 

If  this  sequence  is  monotonic  a  Neville  table  can  be 


. 
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constructed  for  it  too.  Given  an  estimate  g'  of  g  one 
can  obtain  a  better  estimate  for  xc  from  the  sequence 

r  ( n/ ( n  +  g ' ) ) .  (4.1.5) 

n  J 

Alternate  ratios  are  often  used  to  smooth  out  oscilla¬ 
tions  for  loose  packed  lattices.  Fisher  (1962)  suggested 
the  sequences 

=  k  { (n+e)  rn  -  (n+e-2)  rn_2>  (4.1.6) 

for  a  choice  of  small  e . 

In  general  the  series  will  have  non  physical  singu¬ 
larities  too;  the  ratio  sequence  will  converge  more 
rapidly  when  they  are  far  from  the  circle  of  convergence. 

For  series  of  alternating  sign  the  dominant  singu¬ 
larity  will  be  on  the  negative  real  axis  and  for  this  and 
for  the  general  non  positive  series  the  ratio  method  will 
not  yield  the  physical  singularity.  On  the  other  hand  a 
pair  of  singularities  close  to  the  real  axis  could  show 
as  a  consistent  limit  while  the  series  is  well  behaved  on 
the  real  axis  and  has  no  physical  singularity,  but  for 
three  dimensional  cases  we  expect  a  physical  singularity 
so  this  problem  occurs  only  in  the  two  dimensional  case. 
Transformations  of  variable,  usually  bilinear,  have 
been  used  to  bring  the  physical  singularity  in  to  be  the 


* 


closest  to  the  origin  (Wortis  (1969)  ,  Guttman,  Thompson 
(1969)),  with  noticeable  success. 
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(b)  Pade  approximants 

The  (N,M)  Pade  approximant  to  a  power  series  F(x) 
is  the  ratio  of  two  polynomials  of  degree  N  and  M,  whose 
N+M+l  coefficients  are  determined  uniquely  by  the  require¬ 
ment  that  the  coefficients  of  the  power  expansion  of  the 
ratio  be  equal  to  those  of  F(x)  up  to  the  N+M  order. 


[N,M] 


P  (x) 
Q(x) 


,  N 

PQ+P1X  +**PNX 



i+qqx  +..  qMx 


(4.1.7) 


First  introduced  to  critical  phenomena  by  Baker  (1961)  it 
has  proved  an  invaluable  tool  for  the  analysis  of  series 
and  especially  for  non  positive  series  where  the  ratio 
method  does  not  apply. 

In  a  series  of  papers  Baker  has  shown  that  for 
certain  classes  of  functions  the  [m,m]  and  [m-1.,m] 
approximants  form  upper  and  lower  bounds  to  F(x)  for 
real  positive  x  and  that  subsequences  of  the  diagonal 
approximants  converge  to  F(z)  everywhere  in  the  complex 
plane  as  M  ->  °° . 

The  zeroes  of  the  denominator  are  calculated  and 
the  singularity  of  F(x)  shows  up  consistently,  usually 
as  the  smallest  positive  zero.  Obviously  the  approxima¬ 
tion  is  better  if  the  function's  singularity  is  a  simple 


1 
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pole;  therefore  it  is  customary  to  calculate  Pade 
approximants  to 

(i)  ^  log  F (z) , 

or  if  an  estimate  of  y  is  known  to 

(ii)  [F(z)]1//y  f  and 

(iii)  [g^  F(Z)]V^+1  . 

We  have  assumed 

F ( z )  =  (zc  -  z)  ^G(z)  (4.1.8) 

with  G(z)  a  non  singular  function  at  z  .  It  follows 

c 

immediately  that 

(i)  —  log  F ( z )  -  -  ~  log  G(z)  .  (4.1.9) 

c 

So  if  G(z)  has  no  zeroes  near  zc  that  pole  should  appear 
with  good  convergence  in  the  Pade  tables. 

(ii)  (F  (z)  )  1//y+e=  (z  -  z)"1_£  (G(z))1/y+£  .  (4.1.10) 


Here  the  convergence  is  better  as  the  zeroes  of  G(z)  will 
not  interfere.  (iii)  with  an  approximate  y  is: 


f 


3F(z)  ' 
3  z  ' 


1/y+l+s 


1/y+l+s  -1-e 

(-yG (z) )  (Zc_Z^  + 


+ 


( 


3Ch 

3z 


1/y+l+e 


-y/y+1 

(zc-  z) 


(4.1.11) 
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Here  the  convergence  will  be  slower  than  in  (4.1.10) 
because  of  the  second  term  but  still  better  than 
(4.1.9) .  When  the  estimates  for  y  are  not  close  to 
the  real  value  one  expects  a  shift  in  the  estimate 
for  zc  from  Pade  approximants  to  (4.1.10)  and  (4.1.11) 
since 

2 

( z  -  z)“1  +  e  g(z)  =  (z  )~1+e  g(z)[l  +  (z-1)  ~  +  0  (-T)  ]. 

U  L/  Zi  4) 

c  c 


(4.1.12) 


So  to  a  first  approximation  the  pole  of  this  junction 
is  at 


z  (1  +  e) 
c 


varying  the  value  of  e  and  plotting  yc  against  c .  We 
have  a  straight  line, for  e  small  enough  (in  practice  up 
to  rather  large  values)/  whose  intercept  with  the  correct 
value  of  zc  will  give  a  very  good  estimate  of  the 
exponent . 

Another  way  of  estimating  the  exponents  is  to 
calculate  the  Pade  approximant  at  z^  to 

( z - z  )  “  lnF(z)  =  y  +  0  (  z  -  z  )  .  (4.1.13) 

c  cl  z  c 


A  method  which  is  less  sensitive  to  the  choice  of  zc 

was  introduced  by  Baker,  Eve,  Gilbert  and  Rushbrooke 

(1968) ;  the  Pade  approximant  at  z  to 

c 


' 


■ 


■ 


. 


In  F ( z ) 
dz 


(4.1.14) 


d_ 

dz 


In 


dF  ( z ) 
dz 


Y+l  Y+l 

gives  an  estimate  for  J -  .  Since  (4.1.14)  is  J -  at 

Y  Y 

all  z  not  far  from  z  ,  this  is  insensitive  to  the  choice 

c 

of  z  .  A  plot  of  the  Pade  approximants  to  (4.1.14)  gives 

a  straight  line  in  a  y  against  z^  plot.  Its  intercept 

with  the  straight  line  obtained  from  (4.1.10)  yields  a 

simultaneous  estimate  of  y  and  z  . 

c 
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4 . 2  The  static  series  and  their  analysis 

Using  the  graphical  expansions  outlined  in  Chapter 
III,  we  have  obtained  the  initial  parallel  susceptibility 
and  fourth  order  fluctuation  high  temperature  series  for 
the  face-centered  cubic  and  triangular  lattices.  This 
choice  of  lattices  was  made  because  the  close  packed 
lattices  are  known  to  yield  better  behaved  high  tempera¬ 
ture  series. 

For  the  Heisenberg  model  it  is  also  desirable  to 
study  the  staggered  susceptibility  on  loose  packed  lattices 
in  order  to  locate  the  ]\>eel  point,  T  .  However  there  is 
no  reason  to  do  these  calculations  on  the  loose  packed 
lattices  since  for  the  XY  model  the  staggered  parallel  sus¬ 
ceptibility  is  equal  to  the  parallel  ferromagnetic  suscep¬ 
tibility  . 

The  series  obtained  are: 
for  the  face  centered  cubic  lattice, 

v°  =  1  +  6K  +  32K2  +  161. 5K3  +  792. 3K4  +  3823. 924K5  + 

An 

+  18262. 124K6  +  86567. 462K7  +  4 0 284 2 . 0 55 6K8  +  ... 

Y2  =  1  +  24K  +  339K2  +  3656K3  +  33176. 25K4  +  268835K;:)  + 

+  2010189 . 817K6  +  13768318 . 21K7  +  ... 

for  the  triangular  lattice, 

y°  -  1  +  3K  +  7K2  +  13.75K3  +  25.25K4  +  41.0125K5  + 
aii 

+  67.6702K6  +  1 0  9,5  0 1 8  6  K 7  +  170.5454K0  +  ... 


' 

- 


‘ 
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Y2  =  1  +  12K  +  79. 5K2  +  379K3  +  1432. 875K4  +  4621. 5K5  + 
13336. 25K6  +  35115. 962K7  +  ... 

(a)  Comparison  of  susceptibility  and  fluctuation 

In  view  of  the  asymptotic  equivalence  of  the  fluc¬ 
tuation  and  reduced  susceptibility  series  we  wish  to 
compare  the  ratios  of  the  coefficients. 

TABLE  4 . 1 

Ratio  of  fluctuation  over  susceptibility  coefficients 


n 

triangular  lattice 

f.c.c.  lattice 

1 

1 

1 

2 

1.071428571 

1 . 03125 

3 

1 .05454545 

1.01541799 

4 

.965346534 

1 .00435441 

5 

.  971959768 

1 . 00149088 

6 

.  985478098 

1 .00111803 

7 

1.011481242 

1 .00089795 

8 

1 . 026660712 

1 . 01355544 

In  the  case  of  the  f.c.c.  lattice  the  ratios  are 


approaching  unity  ever  more  closely  as  n  increases 
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(except  for  the  last  term) .  The  last  term  breaks  the 
trend  both  here  and  in  the  ratio  plot  and  we  suspect  a 
numerical  error  in  the  susceptibility  series  might  be 
the  cause. 

In  the  triangular  lattice  case  the  ratios  of 
coefficients  of  susceptibility  over  fluctuation  are  more 
erratic.  From  the  result  of  Falk  and  Bruch  y  Y 
at  all  temperatures  therefore  the  ratios  should  be  greater 
or  equal  to  unity  as  n  increases  and  this  seems  to  be  the 
trend . 

(b)  Analysis  of  the  f.c.c.  susceptibility  series 

Figure  4.1  is  the  ratio  plot  for  the  f.c.c.  sus¬ 
ceptibility  series.  We  see  that  the  ratio  plot  is  linear 
(except  for  the  last  term)  and  we  can  estimate 

K_1  =  0.2206  ±  .001  y  =  1.33  ±  0.03  . 

c 

The  above  results  agree  well  with  the  ratio  analysis  of 

the  fluctuation  (Betts,  Elliott  and  Lee  (1970))  who  find 

K  ^  =  0.2206  ±  .0010  and  y  =  1.335  ±  0.02. 
c 

We  apply  the  modified  ratio  method  (Section  4.1) 
to  the  f.c.c.  susceptibility  series. 


r 


, 


' 

■ 
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TABLE  4 . 2 


Ratios  of  f.c.c.  susceptibility  with  Kc=0.24,  y'=1.33 


n 

r 

n 

n (K ' r  -1 )  +1  n 

c  n 

a  /  (n+Y ' -1 ) a 

W  '  n-1 

1 

6 

1.32600 

4 . 51128 

2 

5.33333 

1.35733 

4 . 57798 

3 

5 . 046875 

1.36058 

4 . 54673 

4 

4 . 905882 

1 .33680 

4 . 56023 

5 

4 .826358 

1 .33313 

4 . 52754 

6 

4.775754 

1 .33265 

4.52678 

7 

4 . 740273 

1.33320 

4 . 52686 

Fig . 

4 . 2  exhibits 

. 

(1) 

The  physical 

pole  of  (x°  (K)  )  1//y 

as  a  function  of  y. 

(2) 

The  estimate 

of  y  from  (K  K  ) 

1  c  9K 

log  x„ 

k=k; 

Thi  s 

t 

against  Kc 

is  given  by 

• 

the  vertical  bars. 

The  bars  shrink  as 

approach  K  . 

(3)  The  arrow  at  the  left  hand  side  shows  the  spread  in 
y  calculated  from  Pade  approximants  to: 


w  log  w  x°(K)/w  log  X°(K) 


R '  i 

. 
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The  approximants  to  (-^7)  4  ^  gave  a  straight 

oK 

line  when  plotting  the  pole  against  y.  The  line  had  the 
same  slope  as  (1)  only  with  larger  spread.  We  also 
calculated  Pade  approximants  to  777  log  x (K)  listed  in 

O  I\ 

3  2  ^ 

Table  4.3,  and  to  („-  x/x)  2  listed  in  Table  4.8.  The 

d  W 

results  are  all  consistent. 

TABLE  4 . 3 

Estimates  of  K  for  the  f.c.c.  lattice  from 

c 

Pade  approximants  to  log  X(|(K)  • 


N 

M  =  2 

M  =  3 

3 

II 

M  =  5 

M  =  6 

r- 

II 

S 

0 

.2255 

.2233 

.  2213 

.  2207 

.2206 

.2230+ 

.  005i 

1 

.  2236 

.2095 

.2203 

.  2206 

.  2207 

2 

.  2283 

.  2105 

.2208 

.  2203 

3 

.2201 

.  2208 

.  2205 

4 

.  2201 

.  2201 

5 

.  2206 

/ 


. 


fli 


79 


TABLE  4 . 4 

Estimates  of  K  for  the  f.c.c.  lattice 

c 

4/5 

from  Pade  approximants  to  (x (K) ) 


N  M=2  M  =  3  M  =  4  M  =  5  M=6  M  =  7 


1 

.20017 

.21911 

. 21887 

2 

. 21895 

.  21886 

.21898 

3 

. 21886 

. 21892 

. 21471 
0 . 005i 

4 

. 21911 

. 21722 

.  21892 

5 

.21955 

.21911 

0 

. 21865 

21910 

. 21973 

.  21882 

21820 

.  21909 

21897 

TABLE  4 . 5 


Estimates  of  K  for  the  f.c.c.  lattice 

c 

3 /4 

from  Pade  approximants  to  (x  (K)) 


N 

M  =  2 

M  =  3 

M  =  4 

M  =  5 

M  =  6 

M  =  7 

M  =  8 

0 

.2206 

.2212 

.2212 

.2211 

.  2210 

.2271 

1 

.  2196 

.2216 

.  2212 

.  2212 

.  2208 

.  2211 

2 

.2215 

.  2212 

.2210 

.2210 

.2212 

3 

.2212 

.  2210 

.  2210 

.2210 

4 

.  2212 

.  2210 

.2210 

5  .2207  .2212 

6  .2211 
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TABLE  4 . 6 

Estimates  of  K  for  the  f.c.c.  lattice 

c 

from  Pade  approximants  to  (x,;(K)  ) 


N 

M  =  2 

00 

II 

s 

M  =  4 

M  =  5 

S 

II 

O'! 

M  =  7 

1 

.  22022 

.  22039 

. 22020 

2 

.22045 

.22019 

.22021 

.  22023 

.22022 

3 

. 22021 

. 22021 

.22018 

.  22021 

4 

.22021 

. 22021 

.22021 

5 

.22033 

.22021 

6 

.22017 

TABLE  4 . 7 


Estimates 

of  y  for 

the  f.c.c 

,  lattice  from 

Pade  approximants 

to  (K-Kc) 

A  log  xn(K) 

K=K 

c 

Pade 

K  =  .2200 

K=.  22  02 

K  =.2204 

K  =.2206 

K  =.2208 

K  =  .  221 

K=.  22  2 

c 

c 

c 

c 

c 

c 

c 

(1,3) 

1.293 

1.31391 

1 .327 

1.356 

(2,2) 

1 .294 

1 .3014 

1 .30845 

1 . 31535 

1.32204 

1 . 329 

1.358 

(2,3) 

1.291 

1.2995 

1.30772 

1 . 31579 

1 .32372 

1 .331 

1 .368 

(3,2) 

1.291 

1.2994 

1.30770 

1 . 31579 

1 .32363 

1.331 

1.365 

i — 1 

1 .291 

1 . 2994 

1 .30764 

1.31521 

1 .32213 

1 .328 

1.349 

(2,4) 

1 . 294 

1.3001 

1.30721 

1.31507 

1 .32371 

1 .333 

1.393 

(4,2) 

1.291 

1.2998 

1 .30666 

1.31497 

1.32371 

1.333 

1.386 

(3,3) 

1 .300 

1 .3037 

1.30879 

1.31547 

1.32371 

1 .334 

1 .408 

(1,5) 

1 .232 

11.2706 

1 .29872 

1 .31495 

1.31862 

1.310 

1.119 

(3,4) 

1 .291 

1 . 2995 

1.30772 

1.31580 

1 .32372 

1.331 

1.368 

(4,3) 

1 .291 

1 .2994 

1.30771 

1 . 31579 

1 .32363 

1 . 331 

1.365 

(2,5) 

1.2  91 

1 . 2995 

1 .30765 

1 . 31521 

1 .32213 

1.328 

1 . 350 

(5,2) 

1 .291 

1 .2993 

1.30731 

1 . 31474 

1.32163 

1.328 

1.351 

( 1  r  6  ) 

1.295 

1.3014 

1 .30828 

1 . 31525 

1 .32230 

1 .329 

' 

81 


TABLE  4 . 8 

Estimates  of  K  for  the  f.c.c.  lattice  from 

c 

g  2  i J 

Pade  approximants  to  [  --g--  x(K)/x(K)j 

t)K  li  II 


N 

M  =  2 

M  =  3 

M  =  4 

M  =  5 

M  =  6 

0 

.2217 

.  2206 

.2206 

.2145+  .  06i 

1 

.  2197 

.  2206 

.2206 

2 

.2193 

.  2207 

.  2197 

3 

.  2206 

.2193 

4 

.2206 

(c) 

Analysis 

of  the  fourth  order  fluctuation  series 

for  the  f 

.c.c.  lattice 

The  ratio  plot  is  given  by  the  solid  line  in 

Fig-  4.3.  The  modified  ratios  using  the  known  value 

of  K  are  tabulated  below  and  plotted  in  Fig.  4.3 
c 

(the  dotted  line) . 


TABLE  4 . 9 

Ratios  of  Y2 (K)  for  f.c.c.  lattice,  K^=.221 


n 

c  n 

2 

14.125 

5 .243 

3 

10.785 

5 . 150 

4 

9.074 

5.022 

5 

8.103 

4 .954 

6 

7.477 

4 . 915 

' 
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7 

We  estimate  then  y^=  4.64  ±  0.10  if  we  ignore  the  last  K 

term . 

The  Pade  approximants  to  the  fourth  order  fluc¬ 
tuation  are  not  well  behaved.  A  few  of  the  approximants 
are  tabulated  below.  There  was  no  point  in  plotting 
the  equivalent  of  Fig.  4.1  as  the  spread  in  estimates 
from  different  Pades  was  large,  and  the  against  y 
plot  from  the  estimates  to  [  Y  ^  ( K )  ]  almost  vertical. 

Tables  4.11  and  4.12  are  Pade  approximants  to 

-2—  log  Y„(KX)  and  (KX-  KX)-^—  log  Y  (KX)  with  K= 

3KX  1  C  9KX  1-0. 5KX 

which  seemed  somewhat  better  behaved.  All  the  above  were 
considered  with  the  estimate  from  the  modified  ratio 
method  in  table  4.9. 

TABLE  4.10 

Estimates  of  K  for  f.c.c.  lattice  from 

c 

Pade  approximants  to  (K)  ]  ^ 


N 

M  =  1 

M  =  2 

M  =  3 

M  =  4 

M  =  5  M  =  6 

1 

.2125 

* 

* 

.2376 

.2113  .1919 

2 

.  2269 

* 

* 

.2249 

.  2317 

3 

.  2423 

.  2302 

.  2231 

.2387 

4 

.1908 

.  2127 

.2281 

5 

.2318 

.1853 

6 

* 

*  The 

physical 

pole  did 

not  appear  at  all 

or  was  off  the 

real  axis. 


, 


■ 
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TABLE  4.11 

Estimates  of  K  for  the  f.c.c.  lattice  from 

c 

1  T 

Pade  approximants  to  - — -7  log  Y  (Kx)  where  K=— — - 

3KX  2  1-^KX 


N  M  =  2  M  =  3  M  =  4 


2  .4011  .3634  .3702 

3  .3159  .3697 

4  * 


N 

1 

2 

3 

4 


TABLE  4.12 

Estimates  of  y ^  f°r  the  f-c.c.  lattice  from 

Pade  approximants  to  | (KX-  Kx)  — —  log  Y„(KX)] 

LKX  °  BKX 

where  K  =  — - - 

1-J$KX 


KX= .356 


1 — 1 

II 

s 

M  =  2 

00 

II 

S 

M  =  4 

M  =  5 

4.5188 

6 . 5715 

4 . 0173 

4 .6461 

4 .6300 

5 . 2479 

4 . 3449 

4 . 2448 

4 .6307 

3 . 9788 

4 . 2378 

4 .3037 

4 .7663 

4.7461 

4.7428 


5 


' 
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TABLE  4.13 

Estimates  of  K  for  the  f.c.c.  lattice 

c 

from  Pade  approximants  to  log  Y ^  (K) 


N 

CN 

II 

£ 

M  =  3 

M  =  4 

M  = 

5 

1 

* 

.  2754 

.1938 

.2272 

2 

.  2694 

.  2440 

.  2474 

3 

.164 

.  2473 

4 

* 

TABLE  4 . 

14 

Estimates  of 

y 2  for  the 

f.c.c.  lattice  from 

Pade  approximants  to 

(K- 

-V  3K  l0g 

Y2(K) 

It 

o 

ll 

• 

N 

M  -  2 

M  =  3 

M  =  4 

M  = 

5 

6.367 

4.151 

5.477 

18.436 

5 .473 

4 . 502 

4.24  0 

4 .437 

2 . 969 

4 . 075 

4 


Fig .  4,1 

The  ratio  plot  of  the  initial  parallel 
susceptibility  for  the  f,c,c.  lattice 


Ratios 
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Fig .  4.2 


Estimate  of  K  and  y  for  the  f.c.c.  lattice 

c 

The  solid  line  is  the  estimate  of  Kc  from 


approximants 

to  [X(K)]1/y. 

The  vertical 

bars  are  estimate  of  y  from 

approximants 

to  (K  kc)  gK  log  X  K-K 

c 

The  arrow  on  the  left  is  estimate  of  y  from 

log  x/ff  log  X  at  various  values  of  K. 
o  ■ K  o  i a 
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.219  .2 .lo  .121  *221 

l< 


.2\8 


C 


Fig .  4.3 


The  ratio  plot  for  the  fourth  order  fluctuation 
of  the  f.c.c.  lattice. 

The  solid  line  is  the  ratio  plot  of  the  series. 

The  dotted  line  is  the  plot  of  the  estimates  of  y 

y  =  ( K  r  -  1 )  n  -  1 . 

1  n  c  n 
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( d )  The  triangular  lattice 

The  evidence  for  a  phase  transition  in  the  XY 

model  on  the  triangular  lattice  is  not  conclusive. 

We  present  our  results  below;  a  discussion  is  presented 

in  the  next  section  of  conclusions. 

In  Fig.  4.4  the  plot  of  the  physical  singularity 

o  1  /y 

from  approximants  to  (y  )  against  y  is  given.  On 
the  same  Fig.  4.4  the  estimate  of  y  from  the  approximants 
to  ^  log  x/pr  log  x  is  shown. 

From  Fig.  4.4  we  could  say  that  if  there  is  a 
phase  transition  it  occurs  at  Kc=  .65  ±  .07. 

Fig.  4.5  is  the  ratio  plot. 

TABLE  4.15 

.  .  .  .  •  2 
Ratios  of  triangular  susceptibility  series  Kc=y 

I 

n  r  n  (K  r  -  1 )  +  1 

n  c  n 


3 

2 

2 .33333 

2.1111 

1.96429 

1 . 9286 

1 .83636 

1.8970 

1.62426 

1.4142 

1 .64929 

1.0000 

1.61817 

1.5515 

7 


- 
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TABLE  4.16 


Ratios  of 


_9_ 
9  K 


X„(K> 


for  the  triangular  lattice, 
2 
3 


n 

r 

n 

M  (K '  r  -  1)  +1 
c  n 

1 

4 .66667 

3.1111 

2 

2 . 94642 

2 . 9286 

3 

2.44847 

2 .8969 

4 

2 . 03032 

2.4142 

5 

1 . 98000 

2 .600 

6 

1.88786 

2.551 

3 

The  Pade  approximants  to  log  x (K)  are  not  converged 
at  all. 


TABLE  4.17 


Estimates  of  K  for  the  triangular  lattice 


from 

Pade  approximants 

t°  w  In 

X„<K) 

N 

i — 1 

li 

£ 

M  =  2 

M  =  3 

M  =  4 

LT) 

1! 

M  =  6 

M  =  7 

0 

* 

.5924 

* 

.4814 

* 

.5717 

1 

.  9524 

.7502 

* 

* 

,  5969 

* 

2 

.  5833 

.  8552 

* 

* 

.6552 

3 

* 

* 

.8871 

.  7888 

4 

* 

.3557 

.5988 

5 

.9998 

* 
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TABLE  4.18 


Estimates  of  y  for  the  triangula 

r\ 

Pade  approximants  to  [(K-Kc)  — 


r  lattice  from 
In  X  (K)]  Kc  =  2/3 


N 

M  =  2 

M  =  3 

£ 

II 

M  -  5 

1 

1.8219 

1 . 6724 

1 . 3187 

1 . 6427 

2 

2 .1888 

6 . 2717 

1.4412 

3 

3 .6626 

2 .0076 

4 

1 .2286 

The  analysis  of  the  fourth  order  fluctuation  is 
not  too  conclusive.  The  series  itself  seems  to  diverge 
at  K  =  1.2  ±  0.2  which  could  be  interpreted  both  as 
K  =  1  which  is  just  the  mathematical  radius  of  conver- 
gence  for  a  positive  series  or  as  a  physical  signularity 
at  the  lower  limit. 


TABLE  4.19 

Ratios  of  the  fourth  order  fluctuation  on 
the  triangular  lattice,  =  2/3 

n  r  n  ( K  r  -  1  )  +  1 

n  c  n 


1 

2 

3 

4 

5 

6 


12 

6 .625 
4 .76730 
3.78067 
3.22533 
2.88570 


6 

7 .8333 
7 . 53460 
7 . 08179 
6.75110 
6.54280 
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The  ratio  plots  are  given  in  Fig.  4.6. 


TABLE  4.20 

Estimates  of  y2  for  the  triangular  lattice  from 


Pade 

approximants  to 

(K-Kc>  W  Xog  Y2(K) 

r- 

VO 

vo 

VO 

* 

II 

o 

N 

M  =  2 

M  =  3 

M  -  4 

2 

8.1894 

5.0272 

5 . 0512 

3 

4 . 1519 

5 . 0500 

4 

5.7267 

The  value  of  K  = 

c 

0.6667  was  obtained 

by  a 

better 

analysis  of  the 

fluctuation  series  (Betts  - 

private 

communication) . 

t 


Fig ,  4.4 


Ratio  plot  for  the  parallel  initial  susceptibility 


on  the  triangular  lattice 


/  nn 


r 


T 


T 


T 


9. 


O 


\> 


26 


Ratios 


Fig.  4.5 


Plot  of  against  K  for  the  triangular  lattice. 

The  horizontal  bars  show  the  spread  of  Pade  approximants 
to  [X(K)]i/tr.  The  vertical  bar  is  the  spread  in  as 

estimated  from  -^log fk 
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a 


r 


/ 


•  3 


4- 


■sl/ 


—1 - r— 

.7  .  8 


Fig .  4.6 

Ratio  plot  of  the  fourth  order  fluctuation  series 
on  the  triangular  lattice. 


The  solid  line  is  the  regular  a  /a  Plot* 

The  dotted  line  is  the  plot  of(|T-])  =(K  r  -l)n  against  1/n. 

1  1  1  1 


t 


N>  VP  b-  ^  ^ 


t?6 


Ratios 
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4 . 3  Conclusion 
(a)  Three  dimensions 

We  obtain  from  the  susceptibility  series  on  the 
f.c.c.  lattice  that 

K  =  0.2206  ±  0.001  y  =  1.315  ±  0.04  , 
c 

which  agrees  with  the  results  obtained  from  the  fluctua¬ 
tion  series  by  Betts,  Elliott  and  Lee  (1970) 

K  =  0.2210  ±  0.0005  y  =  1.35  ±  0.03  . 
c 

The  analysis  of  the  fourth  order  fluctuation  yields 
(Ditzian  and  Betts  (1970)) 

y 2  ~  4.64  ±  0.1 

Possible  simple  fractions  consistent  with  our  estimates 
are  : 

4  14 

Y  =  3  Y2  =  X  • 

Therefore  the  gap  parameter  A  is  : 


Using  scaling  relations  (Section  2.4)  we  obtain 
all  other  indices  listed  in  Table  4.21. 


■ 

■ 
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TABLE  4.21 


Index 

Scaling  relation 

Prediction 

Experiments 
for  He 

Y 

4/3 

A 

5/3 

a 

< 

CM 

1 

>- 

+ 

CM 

0 

0 .00010. 003 

3 

A  -  y 

1/3 

0.33310.010 

6 

A/ (A  -  y) 

5 

V 

(  2 A  -  y)/3 

2/3 

0.67  1  0.04 

n 

( 4 A-5y ) / ( 2 A-y ) 

0 

The  experimental  results  are  taken  from  the  following 
experiments : 

The  specific  heat  exponent  a  has  been  measured  by 
Ahlers  (1967) ,  who  measured  Cp  to  great  accuracy 
taking  into  account  gravitational  field  effects  and 
vaporization.  The  temperature  dependence  of  the  super¬ 
fluid  density  in  He4  near  T^  has  been  measured  by  Clow 
and  Reppy  (1966)  .  The  quantity  measured  is  the  angular 
momentum  caused  by  a  persistent  current.  The  result  is 
a  value  for  3  the  exponent  of  the  order  parameter. 

The  value  of  v  has  been  measured  by  Henkel,  Smith 
and  Reppy  (1969)  .  The  experiment  consisted  in  measur¬ 
ing  the  angular  momentum  caused  by  a  persistent  current 
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in  a  film  of  known  thickness  D.  The  angular  momentum 

is  linear  in  D  and  becomes  zero  at  D  .  The  value  of 

o 

D  -D  solid  is  the  healinq  lenqth  E  under  the  assumption 
o  s 

that  near  the  critical  point  only  one  coherence  length 
exists 

E,  cc  (t  ~  T  )“V 
s  c 

We  see  an  excellent  agreement  of  the  XY  model  with 
experiments  on  He^ . 

Domb  (19  )  and  Domb  and  Hunter  (1965)  advanced 

on  Taylor  series  approach  which  predicted  that  the 
exponent  6  be  an  odd  integer.  Our  result  of  6=  5 
agrees  and  also  supports  the  conjecture  that  6=5 
for  all  three  dimensional  models. 

(b)  Two  dimensions 

The  Ising  model  in  two  dimensions  has  not  only  a 
phase  transition  but  the  only  analytically  demonstrated 
phase  transition.  The  other  models  in  two  dimensions 
are  in  a  much  more  dubious  state. 

Mermin  and  Wagner  (1966)  proved  that  spontaneous 
magnetization  cannot  occur  in  one  and  two  dimensions 
for  the  isotropic  Heisenberg  model.  They  show, using 
the  Bogoliubov  inequality  (1962) ,  that  in  two  dimensions 
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|M(H,T) |  <  B  T  In | H | | ^ 
where  B  is  a  constant. 

This  proof,  as  we  mentioned  in  3.3B  does  not 
apply  to  the  XY  model. 

Stanley  and  Kaplan  (1966,  1967)  have  shown  that 
zero  magnetization  at  all  T  >  0  does  not  imply  to  phase 
transition.  The  correlations  can  be  of  a  qualitative 
longer  range  below  some  Tc  and  the  susceptibility  may 
diverge  even  if  M  =  0. 

Stanley  and  Kaplan  (1966,  1967) ,  Moore  (1969) 
show  for  the  classical  Heisenberg  model  and  plane  model 
clear  evidence  of  phase  transition.  Stanley  (1969)  has 
evidence  of  phase  transition  in  two  dimensions  for  the 
quantum  Heisenberg  model. 

The  suggested  behaviour,  Dys  n  (unpublished) , 
Stanley  (1967,  1969),  if 

—  At 

<S  S  >  a  R 
o  R 

such  that  M  =  0  yet  X  =  00  • 

Our  evidence  is  not  really  conclusive.  It  seems 

to  indicate  a  divergence  in  the  susceptibility  at 

K  ~  0.67  ±  0.1  with  an  exponent  y  =  1.5  ±  0.2  which 
c 

agree  with  more  accurate  results  of  the  fluctuation 
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(Betts  -  private  communication)  which  seems  to  diverge 
at  K  =  .67  ±  .01  with  exponent  y  =  1.5  ±  0.02. 


c 


The  fourth  order  fluctuation  series  by  itself 
does  not  help  in  deciding  whether  this  is  a  phase 
transition  or  not  as  it  s  critical  point  is  not  well 
defined  and  it  could  be  that  the  series  converges  up 
to  K  =  1 .  If  there  is  a  phase  transition  at  Kc=  0.667 
then  from  the  modified  ratio  method  y 5.4  ±  0.5. 
Possible  ratios  within  the  error  are 


12 


42 


_  A  15 

then  A  =  -5— 

8  y 


6 
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v  =  8 
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n  =3 
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CHAPTER  V 

DYNAMICAL  CRITICAL  PHENOMENA 

5 . 1  Introduction 

The  critical  region  is  usually  characterized  by 
long  range  and  long  time  fluctuations.  Reaching  equili¬ 
brium  therefore  takes  a  long  time  and  it  is  natural  to 
extend  critical  phenomena  theories  to  cover  dynamic 
properties . 

In  order  to  measure  dynamical  properties  of  the 
system  one  has  to  use  probes  such  as  neutron  scattering 
or  frequency  dependent  magnetic  fields.  While  neutron 
scattering  experiments  depend  on  both  transfer  of  energy 
ftoo  and  transfer  of  momentum  fik ,  frequency  dependent  fields 
have  infinite  wavelength  ,  (k=0) .  For  this  reason 

frequency  dependent  measurement  yields  only  results  that 
do  not  probe  into  the  spatial  dependence  of  the  correla¬ 
tions.  Since  our  calculation  did  not  include  neutron 
scattering  we  will  not  discuss  it. 

In  the  Ising  model  correlations  of  the  z  components 
of  the  spins  can  only  have  space  dependence  as  each  cv 
commutes  with  the  Hamiltonian  and  therefore  is  constant 
in  time.  Fisher  and  Burford  (1967))  calculated  the  zz 
correlations ,  which  give  the  elastic  scattering  of  neutrons 
from  an  Ising  model. 


- 

. 


. 


. 
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Allan  and  Betts  (1968)  using  the  Kubo  linear 
response  method  calculated  exactly  the  frequency  depen- 

XX 

dent  susceptibility  x  (w)  anc^  inelastic  neutron  scattering 
from  a  spin  Ising  model  in  2  dimensions. 

Essam  and  Garelick  (1968)  obtained  the  exact  solution 
for  all  spins  and  in  presence  of  parallel  fields,  their 
solution  was  based  on  solving  the  Green's  function  hierarchy 
of  operations  of  motions  without  need  for  approximations 
other  than  linear  response. 

Allan  and  Betts  (1968)  noticed  that  the  limit  of 
X  ( co )  as  03  goes  to  zero  was  equal  to  xT  (0)  in  the  case 
of  the  honeycomb  lattice  but  not  the  square  lattice. 

The  Kubo  "adiabatic",  or  quasistatic  susceptibility  is 
not  equal  to  the  isothermal  susceptibility  on  the  square 
lattice  because  the  system  is  not  ergodic  in  the  absence 
of  a  field  in  the  x  direction. 

As  shown  by  Essm  and  Garelick  (1968)  M  =  -  £  P  E*  (H) 

r 

where  E^.  (H)  is  the  first  derivative  of  the  r 1  th  energy 
level  with  respect  to  the  field  and  P  (T) ,  the  probabi¬ 
lity  of  finding  a  number  of  the  ensemble  in  state  |r> 

Xxx  =  ^  =  _E  P  e;(hx)  .-  e  ^  e;(hx) 

T  3  HX  T  r  d  T 

3  P  , 

=  X  (w=0)  -  Z  (^jp)  Er(HX) 


(5.1.1) 


' 


■ 
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and  e' (Hx=  0)  is  not  zero  when  there  is  degeneracy,  that  is 
r 

different  spin  configurations  with  same  energy. 

Allan  and  Betts  (1968)  pointed  out  that  spins  which 
have  an  equal  number  of  neighbours  pointing  up  as  down, 
are  essentially  free.  This  causes  the  degeneracy 
mentioned  above.  We  want  to  note  here  that  for  the  XY 
model  in  any  lattice  even  1  dimensional  there  is  no  such 
degeneracy  when  we  deal  with  the  magnetization  in  the  x 
direction,  which  is  parallel  to  the  interaction.  As 
discussed  in  the  next  section  the  magnetization  in  the 
z  direction  is  non-ergodic  at  least  in  the  1-dimensional 
case . 

For  the  XY  model  we  had  to  use  series  expansion 
methods  for  the  time,  or  frequency,  variable  as  well 
as  for  the  temperature  while  the  Ising  model  results 
mentioned  above  had  closed  formulae  for  the  time  variable. 
Series  expansions  were  only  needed  for  the  static  spin 
correlations  of  the  spin  at  the  origin  with  k  of  its 
nearest  neighbours,  k  ^  q. 

Another  subject  we  want  to  mention  is  the  critical 
slowing  down  for  dynamical  susceptibility.  Suzuki  (1969 
a,b)  calculated  the  dynamical  susceptibility  for  the 
kinetic  Ising  model  in  two  dimensions (and  for  the  XY 
chain  which  we  discuss  in  the  next  section) .  We  quote 


/ 
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here  some  of  his  results  as  this  was  the  only  calcula¬ 
tion  of  a  dynamic  susceptibility  for  a  physically 
plausible  model  before  our  work,  and  we  use  his  defini¬ 
tions  and  terminology. 

The  kinetic  Ising  model  has  transition  probabilities 
defined  by 

oj  .  (a  )  =  ~  (1  -  a  .  tgh  3  Z  J.,  a,)  .  (5.1.2) 

3  o  /t  j  ^  3K  K 

t  indicates  the  relaxation  time  of  a  free  spin  interacting 
with  a  heat  bath.  This  system  has  a  Louville  operator 

—  Z  w^(q^)  ^ ~  ^ k ^ 
k 

where  P,  flips  the  k'th  spin.  The  dynamical  susceptibility 
can  be  expanded  in  a  moment  series: 

x  (u)  =  N  U  i  { J—  *  - J  +  ...}  (5.1.3) 

1UT  (iut) 

where 

u  =  — <M  M>  •  (5.1.4) 

Nm 

All  moments  u  are  positive  and  remain  finite  at  the 

n 

critical  point.  Suzuki  proved 


< 


r 
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Nm2$jj 
(cot  )  2  + 


2 


-<  Re  x  (co ) 


Nm23li  2 
(cot  )  2 


(5.1.5) 


and  therefore  postulates  that  the  susceptibility  behaves 
like 


X(to)  ~  - - - — -  (5.1.6) 

ioo  +  sY  f  (—  ,  e ) 
eY 

and  f(0,e)  can  be  singular  at  e=0 .  Expanding  in  powers 
of  co 


n 


X  (co )  =  N$  Z  (-icox  ) 


and  calculation  gave 


-a. 


a  (e )  ~  e 
n 


n 


anM  +  n  As 


(5.1.7) 


(5.1.8) 


where 


A 


s 


(5.1.9) 


and  A  was  proved  by  Abe  and  Hatano  (1969) . 

s 

of  critical  slowing  down  A  is  defined  from 

s 


/ 

o 


<M (t)M> 

<m2> 


dt 


The  index 


e 


(5.1.10) 


' 

. 
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and  is  found  to  be  the  same  as  in  (5.1.9) .  Suzuki 
defines  the  index  of  critical  slowing  down  of  the 
self  correlation  (see  5.4) 

°o  -A 

/  <c,  (t)  a-,  >  dt  ~  e  .  (5.1.11) 

0  1  1 

Classical  theories  predict  A^  =  y.  Suzuki  defines 
^(e)  by 

tM  =  x(0)  RM1(e)  (5.1.12) 

and  states  a  proposed  law  of  similarity  as  follows. 

The  relaxation  time  singularity  consists  of  the  direct 
or  thermodynamic  critical  slowing  down  due  to  the 
susceptibility  and  the  indirect, due  to  R(e)/  called 
induced  critical  slow  down  (or  speed  up  as  the  case 
may  be) .  Relaxation  times  can  be  defined  for  various 
quantities  not  only  the  magnetization.  The  similarity 
law  states  that  the  indirect  critical  slow  down  is  the 
same  for  some  critical  variables .  Variables  for  which 
it  holds  are  called  similar  variables. 

For  the  kinetic  Ising  model  in  2  dimensions  Suzuki 
expects  that  the  magnetization  and  the  energy  are  similar 


variables . 


' 


■ 
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T 


E 


tem 


<5E  X 

<6E  i 
<  M  7 


6E  >/  <  (6E) 2> 
M  >/ <  M6E>  ~ 
M  >/<  M2>  - 


(5.1.13) 

(5.1.14) 


(5.1.15) 


and  we  know  <(SE)2>  ~  z  a  and  <M6E>  ~  |e|^  \  Only  A^ 

is  calculated  directly  but  Suzuki  using  the  relation 


< (AM  +  6E)  i  (AM  +  6E) >  >>  0  (5.1.16) 

O Ls 

proves  that 

2  AME  ”  (AM+AE^^a  +  2  3  +  Y  “  2  (5.1.17) 

for  real  X.  Similarity  and  (5.1.17)  yield  a  +  23  +  Y  >-  2 

which  is  valid  assuming  scaling, Suzuki  obtains  that 

Am=2.00  ±  0.05.  Similarity  claims  A£  -  a  =  AM  -  y  and 

7  i 

since  a  =  0  y  =  j  this  yields  A£  =  This  is  consistent 

with  computer  simulations. 

We  calculated  a  series  for  x  (w)  by  methods 

explained  in  sections  (5.4)  and  (5.5)  and  estimated  the 

critical  index  of  slowing  down  A  obtaining  that  in  this 

s 

case  it  is  actually  a  speeding  up,  A 

o 


<  y. 


' 
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and  gH  (°°,  0)  =  g  (°°,  0)  .  y  =  d+x  where  d  is  the  dimen¬ 
sionality.  The  indices  n  and  v  were  defined  in  (2.1.8) 
and  (2.1.9)  respectively.  The  Fourier  transform  of 
(2.1.9)  is 

F (k ,  e=0 )  «  ky  Z 

r  (k=0 ,  c )  cc  e  ^  oc 

as  E,  ~  e"1//v  .  From  (5.2.4)  we  see 

T(k,  £)  a  E, ~Y  [1  +  ...]  for  k^  <<  1  (5.2.7) 

if  r  (k=0 ,  e)  is  finite  for  ~  /  0.  Similarly  if  F(k,  e=0) 
is  finite  for  k  ^  0 

T(k,  O  «  ky  [1  +  ...]  k£  >>  1  .  (5.2.8) 

Comparing  (5.2.5),  (5.2.6)  with  (5.2.7)  and  (5.2.8),  the 

scaling  law 

y  =  —  2  +  r]  =  —  y/v  (5.2.9) 

results.  Adding  assumptions  on  the  structure  of  r(k,£,H), 
all  scaling  laws  (2.4.24)  predicted  by  the  Kadanoff  cons¬ 
truction  are  obtained.  Returning  to  r(r,£),  since 


(5.2.5) 

(5.2.6) 
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x  =  y  -  d 

T(r,  £)  =  rn  2  d  g(r/£)  (5.2.10) 


The  dynamical  scaling  assumption  is  a  generaliza¬ 
tion  of  (5.2.4)  to  the  frequency  dependent  correlation. 
r^(k,  go)  is  defined  by 

d^k  ,  i(k.r-Got) 

rr  (r ,  t)  =  /  - - rr(k,u>)  (5.2.11) 

S  -  (2,0 3  2*  5  - 

where 


I^(r,t)  =  <{S(r,t)  -  <S(r,t)>  ,  S(0,0)  -  <S(0,0)>}>  . 

(5.2.12) 


The  curly  brackets  are  anticommutator.  B, ,  the  correla¬ 
tion  length,  gives  the  dependenos  on  the  temperature. 

The  scaling  assumption  is 


(k  ,go) 


2 77  f  (k) 


GO 

fl(kg) 


) 


kz  n(kC) 


(5.2.13) 


where 


r?(k) 


oo 

/ 

— oo 


doo 
2  7T 


(k , oo) 


(5.2.14) 


and  (5.2.14)  implies  then 


—  00 


f  f  (x)  dx  =  1  . 


(5.2.15) 
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ft  is  determined  from  the  constraint 


1 

/  f,  (x)  dx 
-1 


1 

2  * 


(5.2.16) 


The  function  kZ  ft(k£)  is  called  the  characteristic 

M 

frequency  co^  (k)  .  The  time  dependent  correlation  is 
scaled  as 


rr (k)  f (t.kz  ft (k£) ) 

rr(k,t)  «  _5 -  .  (5.2.17) 

4  kz  ft( k£) 


For  the  isotropic  Heisenberg  model  the  Landau- 
Lifshitz  (1935)  formula  gives 

(Is ^  =  ^ 

where 


<  M  > 


and  p  is  a  stiffness  constant  proportional  to  £ 
s 

From  this  Halperin  and  Hohenberg  deduce  to  first  order 
z  =  3  -  3/v 

z  cannot  be  determined  from  thermodynamic  considerations 
in  the  anisotropic  case. 


' 

•  n  :,¥  ' !  1  1 


■ 
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We  calculated  <a^(t)  a^>  ,  the  self  correlation. 
The  NMR  line  width  AfflR  is  given  by  (Heller  (1966)): 


anmr  *  1  <  al(t)  °1  >  dt 


(5.2.18) 


00 


/  r  (r=0 ,  t)  dt 
0  ^ 


oc  /  dJk  T  (k,u)  =0) 
0  ^ 


oo  k  dk  rr  (k)  f.  r  (0) 

cc  /  - i ^ - 


0 


k  fi(k£) 


/  dk  k 
0 


2-z+n-2  g(k?)  fk^(0) 


n(k?) 


cc  £z  n  1  /  dx  F  (x) 


0 


Since  E  ~  e 


-v 


A  <x£ 
NMR 


-v  (z-n-l) 


(5.2.19) 


The  assumption  that  (5.1.13)  holds  for  the  order  para¬ 
meter  is  restricted  scaling.  Extended  scaling  is  when 
(5.2.13)  holds  for  other  operators  also. 

A 

3 

Assuming  that  the  frequency  scales  like  co/e 
yields,  using  (5.2.13),  that 

z  =  v"1  A  .  (5.2.20) 
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5 . 3  Dynamical  results  for  the  one  dimensional  XY  model 

Many  authors  have  studied  the  time  dependent  corre¬ 
lations  in  the  one  dimensional  XY  model.  (a)  Niemejer 
(1967,  1968),  (b)  Katsura,  Horiguchi  and  Suzuki  (1969), 

(c)  Suzuki  (1969),  and  (d)  Barouch,  McCoy  and  Dresden 
(1970)  . 

(a)  and  (d)  were  concerned  with  the  exact  solution 
of  the  anisotropic  XY  model  which  can  be  obtained  for 
arbitrary  fields  in  the  z  direction.  Of  course  specific 
cases  and  asymptotic  limits  had  to  be  taken  to  see  what 
those  solutions  mean.  (b)  obtained  linear  response  for 
the  isotropic  XY  model  by  the  two  time  Green's  function 
method.  (c)  discussed  the  results  for  the  anisotropic 
XY  model  terms  of  critical  slowing  down. 

In  (a)  the  static  result  obtained  was  that  the 
one  dimensional  XY  model  has  a  divergence  at  T=0,  h=h  , 


^\2> 


He 

V7Z(h'  T=  o) 


Kc 

<  T=  o  )> 


- 

. 


' 
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This  makes  the  order  in  which  limits  are  taken  to  obtain 
asymptotic  behaviour  crucial  .  The  susceptibility  behaves 
like 


XT=0  ~  lo9  ( I h-1 1  ) 


where 


h 


(5.3.1) 


y  being  the  anisotropy  parameter.  The  points  of  interest 
to  us  were  in  (d)  and  (c) . 

In  (d)  the  authors  found  that  the  asymptotic 
behaviour  of  M(t)  is  the  same  for  large  times  for  all 
H,only  the  coefficients  depending  on  the  details.  There 
is  no  approach  to  equilibrium  no  matter  how  slowly  the 
field  varies.  The  easiest- example  they  provide  is 

a  t  0 

H(t)  =0  t  >  0  . 

V. 

The  plot  of  M  (t=~  ,  T=0 )  as  a  function  of  a  looks  like 

z 


with  M  definitely  saturating  at  a  value  less  than  Jg, 
z 


i 

. 

. 
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where  \  is  the  saturation  value  of  the  initial  magneti¬ 
zation  . 

Suzuki  obtains  that  the  dynamical  susceptibility 
behaves  like 


Re  h  >  T=0)  ~  -log  |  a)  |  (5.3.2) 

c 

and  despite  lack  of  ergodicity 
lim  (M,  M ( t ) )  =  <M> 2 

t+  oo 

so  he  can  define  the  relaxation  time  and  obtains 


2  “  (SM,  6M  (t)  )  t  ..  u  2AM 

tm  =  - (mTm) - dt  h 


where  A„  =  1 . 


For  the  partial  energy  defined  by 


E  =  4  E  S*  S* 

l  i+l 
i 


(5.3.3) 


(5.3.4) 


he  obtains  A  =  1  and  since  y  =  0  similarity  holds  for 

ill  lli 

this  partial  energy  and  the  susceptibility. 

We  note  that  here  A^  >  y^but  here  because  of  the 
singular  character  of  the  model  the  temperature  was 
kept  fixed  at  T=0  and  the  magnetic  field  is  the  variable^ 
which  is  quite  different  from  our  procedure.  Also  these 
results  were  obtained  by  Suzuki  for  the  anisotropic  XY 
model  as  Re  x  (k=0 ,  w)  =  0, according  to  his  calculations, 
when  the  anisotropy  vanishes. 


' 
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5 . 4  The  Kubo  linear  response  theory 


The  Kubo  formalism  (Kubo  1957)  is  a  general  approach 
to  the  problem  of  time  dependent  disturbances.  The  results 
are  perturbation  series  in  the  disturbance.  We  shall,  as 
most  users  of  the  formalism  do,  keep  only  the  term  linear 
in  the  disturbance. 

The  system  is  assumed  to  be  in  equilibrium  without 
disturbances  at  t  =  -°°.  Its  properties  are  then  described 
by  a  density  matrix  p  .  For  a  canonical  ensemble 


P 


o 


Tr  e 


3^o 


(5.4.1) 


A  disturbance  is  introduced;  the  external  force  F(t)  is 
conjugate  to  the  observable  A.  The  perturbed  Hamiltonian 
is  now 

U  =  M0  +  it'  =K0  -  A.F(t)  .  (5.4.2) 

The  time  evolution  of  the  system  is  given  by  the 
equation  of  motion  of  the  density  matrix 


cl  p(t) 


_1_ 

ifi 


[K  ,  p  (t)  ] 


(5.4.3) 


The  solution  of  (5.4.3)  after  linearization  is: 


i  (t-t 1  )XQ  i(t-t');£0 


p(t)  =  pQ  +  ih  /  e 

—  00 


[^ '  P0^l 


h 


dt  . 


(5.4.4) 


■ 


. 
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The  response  of  an  observable  B  to  the  disturbance  is 


AB  =  Tr  (p(t)  -  pQ)  B  . 


AB  is  the  macroscopic  difference  of  the  average  of  B  at 
time  t  from  its  equilibrium  average.  It  follows  that: 


Tr 


A]B  (t-t '  )F  ( t '  )dt ' 


(5.4.5) 


where 


B  ( t) 


B 


The  linear  response  function  <f>  ,  which  is  the 

response  of  the  observable  B  to  a  disturbance  conjugate 
to  A  is  defined  by 


AB  (t) 


t 

/ 


$B^ ( t-t 1 ) F ( t 1 )  dt' 


(5.4.6) 


It  follows  then,  that 

«i>BA<t)  =  nr  Tr  p0  tA'  B(td  •  (5.4.7) 


Many  disturbances  of  interest  can  be  expanded  in 
Fourier  series.  We  therefore  seek  the  response  to  a 
harmonic  force.  It  is  convenient  to  take  F  a s 


F  (t ) 


Re  F 


iwt+et 


(5.4.8) 
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where  e  is  a  small  positive  factor  which  ensures  that 
F  ( — oo)  =0  as  an  initial  unperturbed  condition.  This  is 
a  slow  quasistatic  switching  on  of  the  force. 

The  frequency  dependent  susceptibility  is  defined 

by 

AB(t)  =  Re  XBA(^)  F0  e1Wt  •  (5.4.9) 

The  result  of  combining  (5.4.9)  and  (5.4.6)  is: 

Xra(w)  =  lim  /  <|>_.(t)  e~la)t  dt  .  (5.4.10) 

BA  6+0+  0  A 

The  response  function  is  real  and  the  susceptibility 
satisfies  the  following  symmetry  relations 

(1)  Fe  XBA(u)  =  Re  XBA(“w) 

(2)  Im  XBA(^)  =  ~Im  XBA(-w) 


XBA^'  CA  £B  ^AB^W'  F^  ' 


where  e  ,e  are  +1  or  -1  as  A  and  B  are  respectively  even 
A  B 

or  odd  with  respect  to  time  reversal.  Relation  (3)  is 
known  as  the  Onsager  relation  (1931)  . 

For  short  times  the  response  function  can  be  expanded 
in  an  asymptotic  power  series .  Using  (3-3.17) 


*BA(t)=4Tr5,°  [a’(b+7%b]+  (ir)2[rf’[w,B]] +- •  •  )]  (5-4-11) 


=  53  f  t 
n 


n 


' 


. 
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This  series  is  either  odd  or  even  in  t  and  in  our  magnetic 
case  (J)^  is  odd.  Substituting  (5.4.11)  in  (5.4.9)  and 
using 


.  m  cx  ioox  *  m  ^m ,  ■  \  . 

/  x  e  e  dx  =  tti  6  (oo-ie)  + 

0 


ml 


(ioo+e ) 


m+1 


(5.4.12) 


we  obtain  for  6^..  odd  in  t 
'  yBA  ' 

m+1 


Re  xBA(aj)  =  E 


m!  (-1) 


oo 


m+1 


m 


(5.4.13) 


Im  xBA(w)  =  ir(i)m  «m(“)  f 


m 


The  initial  values  of  the  response  function  (since 

f  =  -A-  <j>  ^  (0))  determine  the  coefficients  of  the  power 
m  ml  Y 

series  in  —  for  the  susceptibility  and  the  moments  for 
oo 

the  dissipative  part.  The  quantity  we  are  calculating 
here  is  yXX(w)f  that  is  the  force  we  apply  is  HXelwt  and 
both  operators  A  and  B  are  m  .  The  response  function  is: 

♦xx(t)  =  4  Tr  po  fV  Mx(t)]  (5.4.14) 

and  the  susceptibility  follows  from  (5.4.11)  and  (5.4.13). 
Obviously  in  this  case  6  (t)  =  -(}>  (~t)  .  It  can  be  noted, 

d>  (t)  =  4>  (t) 

Txx  Yyy 


that 
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and 

<p  (t)  =  d>  (t)  =  (f>  (t)  =  cj)  (t) 

XY  Yyx  Yxz  Yyz 

=  <p  (t)  =  <j>  (t)  =  0  . 

But  for  neutron  scattering  where  one  needs  the  x ( ^ r  w ) 
the  zz  part  does  not  vanish. 

We  want  to  obtain  the  zero  frequency  susceptibility. 
The  identity 

“  _  p  v  3 

[A,  e  °]  =  V-  e  ^  f  A(-iRX)  dA  (5.4.15) 

1  0 

was  proved  by  Kubo  (1957)  and  independently  by  Feynman 
(1948) •  Using  (5.4.15)  and  integrating  by  parts  we 
obtain  the  following  expression  in  the  zero  frequency 
limit 

B  o  o 

XBA(0)  =  f  Tr  po  A("ifiX)B  dA  -  3  Tr  pQ  A°BU  . 

(5.4.16) 


XD7V(0)  is  the  quasistatic  (also  called  adiabatic)  suscep- 
dA 

tibility  and 

1  T 

A°  =  1 im  —  /  A ( t )  dt  . 

0 


is 


The  isothermal  susceptibility  on  the  other  hand 


Tr  p  B  -  Tr  p  B 
^  o 


F 


(5.4.17) 


Idt 


' 
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which  can  be  shown  (Kubo  1957)  to  be  equal  to 


X>> 


3 

f  Tr  p  A(-ifU)B  dA 
0 


-  B<A><  B> 


(5.4.18) 


in  the  linear  approximation. 

The  second  term  in  (5.4.16)  was  obtained  from 

3 

lim  /  Tr  p  A  (-ill A  )  B  ( t)  dA 
t-^°°  0 


by  noting  that  if  this  limit  exists  it  is  equal  to 


lim  ~  /  dt  /  Tr  p  A(-ihA)  B(t)  dA  - 
t  o 

t->°°  0  0 

3 

=  /  Tr  p  A(-ifiA)  B  dA  . 

A  *0  o 


This  limiting  process  picks  up  the  zero  frequency  compo¬ 
nent  of  B  f  Bq ,  which  is  diagonal  with  respect  to  <ftQ . 

Now 


Tr  p  Q  B  =  Z  plj  Qjk  Bkl  = 

Ho  o  .  .  Ko  o 

1,] 

=  £  p11  Q11  B11  =  Tr  p  0  B 
.  Ko  o  o  o  o 

l 

and  the  diagonal  part  of  A(-ifiA)  is  just  A°  since 
(e  A  e  )±i  =  Au. 

For  the  isothermal  and  adiabatic  susceptibility  to 


be  equal  it  is  sufficient  that 

lim  <AB  (t ) >  =  <  A>< B> 
t->°° 


(5.4.19) 


A  system  which  satisfies  (5.4.19)  is  called  ergodic.  For 
a  discussion  of  the  ergodicity  of  our  system,  see  Section 


(5.1  ) . 


. 


' 


- 
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5 . 5  The  high  temperature  series  expansion  for  the 
frequency  dependent  susceptibility 

As  we  saw  in  the  previous  section,  the  susceptibility 
can  be  expanded  in  an  asymptotic  series  in  powers  of  co  ^  . 
The  coefficients  of  this  expansion  are  determined  by  the 
time  derivatives  of  the  response  function  at  zero  time. 
Since 
t 

<M  (t)>  =  /  <f>  (t-t1)  H  (t1)  dt '  (5.5.1) 

x  1  rxx  x 


we  expect  <j>  to  be  linear  in  N,  the  number  of  lattice 

XX 

sites.  Therefore  in  the  terminology  of  Sec.  3.2 
4>xx  =  ’  {linear  part  in  N  of  -jr<  [Mx  ,Mx  (t)]  >} 

=  -i  Tr  eKP  [Mx,  Mx(t)]  .  (5.5.2) 


Defining  oj  =  J/h  we  have 

-iPto  t  iPw  t 

d>  (t)  =  -i  Tr*  eK  Qm  ,  e  M  e  ]  .  (5.5.3) 

Yxx  ft  L  x  x  J 


We  expand  all  exponentials  in  (5.5.3)  in  power  series 
of  P  and  have  a  double  series  in  K  and  t 


<f> 


xx 


.  v2 

lm 

4ft 


Tr1 


OO 

z 


m, n , q=0 


(iQ)0t)m+q 
n  !  ml  q I 


I  (Pna* (-P) moxPq-  Pn(-P)mc^pV) 

i»  j  1  3 


X 


t 


(5.5.4) 
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where  m'  is  the  magnetic  moment  of  the  atoms. 

We  rearrange  the  sum,  use  the  cyclic  property  of 
the  trace  and  interchange  m  and  q  in  the  first  term  of 
(5.5.4)  to  obtain: 


4> 


xx 


Tr* 


E 

n,  s 


(iwQt) 


x 


s 

E 

m=0 


0 

m 


E 

if  j 


,n+m  x  _  s -m  x  ,  ,  %  s  m  .  -1  \  m*  /  c  c  c  \ 

°i  P  aj  “  (-D  )*  (5.5.5) 


For  even  s  [(-1) 


s-m 


(-l)m]  vanishes.  Therefore, 


4> 


.  *2 
lm 

xx  2h 


Tr5 


L  nl 
n 


odd 


( io3Q  t ) 
s  I 


X  E  (s)  (-l)m  E  Pn+m  oX  PS“m  0X 
m=0  m  i,j  1  D 


(5.5.6) 


x  f 

Substituting  ai  =  ai  +  ai 


.  v2  oo  n 

<t>  (t)  =  E  — 7-  E  , 

xx  2h  n=0  n*  s  odd 


(iu)Qt) 
s '. 


x 


X  X  (s)  (-l)m  Tr*  [  X  (Pn+m  at  Ps'ra  a,  +  Pn+m  a.  P5’”  at] 
m=0  m  i,j  13  13 


(5.5.7) 


We  have  in  (5.5.7)  an  expression  similar  to  the  one  we  had 
for  the  static  susceptibility  (3.3.10)  with  an  added  com¬ 
plication.  Now  a  different  factor  multiplies  different 


123 


vertical  orderings  of  the  arrows. 

Instead  of  one  vertical  weight  we  now  calculate 
£+1  weights,  where  Z  is  the  number  of  solid  arrows  in 
the  graph.  As  before  the  tip  of  the  dotted  i->j  arrow 
is  fixed  -  at  the  zero  level.  The  tail  is  fixed  in 
sequence  at  each  of  the  levels  k  and  the  number  of 
allowed  orderings  of  the  other  arrows  is  v(k-2) .  For 
a  fixed  n  and  s 

Z  (s) (-l)m  Tr*  Pn+m  a+  Ps"m  a.  = 
n  m  1  1 

m=0  J 


I 

i 

^s+n+1 


l  (S)  ,  j-m  (gs+n>L)h(g')v(s-m) 
m=0  m  2V 


(5.5.8) 


where  v  is  the  number  of  vertices. 

As  in  the  susceptibility  the  graphs  were  either 
"bubbles"  for  i=  j ,  or  fluctuation  type  graphs  with  a 
dotted  arrow  i-j  unrestricted  in  length.  It  can  be 
noted  though  that  when  v(k)=v  for  all  k,  we  have  in 
(5.5.9) 

Z  O  (-Dm  v  =  0 

m 


so  graphs  where  the  i-j  arrow  has  a  vertex  of  order  2 
do  not  contribute,  which  reduces  greatly  the  number  of 


■ 
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fluctuation  type  graphs.  There  is  no  need  to  calculate 
all  v(k),  since  using  the  cyclic  property  of  the  trace 
we  can  obtain  that 

v(k)  =  v ( Z  -  k)  .  (5.5.9) 

The  equivalent  of  Theorem  I  is  a  bit  more  complicated 

here . 

Theorem  la 

Let  g  be  a  disconnected  graph  of  Z  solid  arrows  and 
one  dotted  arrow.  There  are  two  subgraphs,  g^  and  g^. 
g^  has  Z-^  solid  arrows  and  a  dotted  one  and  g^  has  J 
solid  arrows.  Then: 

v  (k)  =  v2  E  v1(i)  (k^1)  (J  ^1)  (5.5.10) 

i=0  1 

where  (i)  is  the  vertical  weight  of  g^  with  the  dotted 
tail  at  the  ith  level.  v(k)  is  the  vertical  weight  for  the 
whole  graph  g  with  tail  of  the  dotted  arrow  at  the  kth 
level.  v^  is  the  vertical  weight  of  as  a  partition 
function  graph. 

Proof 

We  keep  the  dotted  tail  fixed  at  the  kth  level. 

Among  the  k-1  levels  between  the  tip  and  tail  of  the 
dotted  arrows  i  will  be  filled  by  arrows  from  g^. 


The 
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number  i  can  run  from  zero  to  k-1  and 

k-1 

v(k)  =  v?  Z  v, (i+1)  X. 

A  i=0 

where  is  the  number  of  ways  £  arrows,  £^  e  and 

£2  e  g2  >  can  be  assigned  to  £  levels  in  such  a  way  as 
to  have  i  arrows  belonging  to  g  assigned  to  the  first 
k-1  levels. 

The  result  (5.5.10)  follows  when  we  note 

„  _  (k-1)!  ( £-k+l ) 

1  il(k-l-i)  (£1-i)l  (£2-k+l-i)!  * 

The  first  factor  is  choosing  i  levels  out  of  k-1  to  fill 
with  g1*  The  second  factor  is  the  number  of  ways  of 
assigning  £-k+l  levels  so  that  £1~i  are  occupied  by  arrows 
of  and  the  rest  by  arrows  of  g2. 

The  equivalent  of  Theorem  II  holds  too. 

Theorem  I la 

Let  a  connected  graph  g  of  £  solid  arrows  and  a 
dotted  arrow  be  composed  of  two  subgraphs,  g1  of  £^  solid 
arrows  and  the  dotted  arrow  and  g2  of  ^  s°lid  arrows 
such  that  g1  and  g2  have  m  common  vertices  of  order  2  in 
g  and  no  other  common  vertices.  Then  the  vertical  weights 
v(k)  of  g  are  given  by  the  same  expression  (5.5.10)  as  in 


Theorem  la. 


* 


Proof 


Exactly  as  the  proof  of  Theorem  II  followed  that 
of  Theorem  I.  We  have  no  restrictions  on  vertices  of 
order  2,  so  they  can  be  cut  to  give  two  separated 
graphs . 

A  numerical  example  to  clarify  the  meaning  of  the 
theorem:  The  graph  g  is  composed  of  g^  and  g2 


1=6,  &i  =  £2  =  v2  , the  vertical  weight  of  g2 

alone/ is  3  1  =6.  v^(i)  are  the  vertical  weights  for  g^ 

where  i-1  levels  separate  head  and  tail  of  the  dotted 
arrow.  Obviously  v-^(l)  =  v-^(4)  =  1  ,  v-^(2)  =  v-^(3)  =  0. 

Let  us  calculate  v(3),  the  vertical  weight  of  g  when  2 
levels  separate  head  and  tail  of  dotted  arrow. 

Substituting  in  (5.5.10) 

v  (3)  =  6.1.  (q)  (3)  =  24  . 

With  some  labour  the  same  result  can  be  obtained  directly. 
Theorem  Ilia 

Let  a  connected  graph  g  of  l  solid  arrows  and  a 
dotted  arrow  be  composed  of  two  subgraphs  g^  and  g2  of 


' 

1 
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and  solid  arrows  respectively  and  a  third  subgraph 
consisting  of  the  single  dotted  arrow.  The  subgraphs 


g^  and  g^  have  in  common  m  vertices  of  order  2  in  g  and 


no  other  common  vertices.  Then  the  vertical  weights  are 
given  by: 


k 


(5.5.11) 


v(k)  =  Z 
i=0 


The  dotted  arrow  belongs  to  both  g  and  g^  for  the  purpose 
of  calculating  the  v^i+1),  v2  (k-i)  .  The  proof  follows  exactly 
the  proof  of  la.  After  assigning  i  levels  from  the  first 
k-1  levels  to  g^ ,  we  assign  the  (k-l-i)  levels  left  to  g2* 

The  dotted  arrow  belongs  to  both,  so  the  vertical  weights 
of  both  are  dependent  on  i. 

In  the  actual  calculation  we  did  not  go  far  enough 
to  need  the  theorems  ytheir  only  use  was  for  checking  the 
computer  programs. 

Using  the  fact  that  vertices  of  order  2  introduce 
no  added  restriction  the  dotted  arrow  was  broken  by  adding 
a  vertex  of  order  2  and  the  two  new  arrows  were  put  at 
the  head  of  the  vector  describing  the  graph. 


In  Appendix  I  the  graph 


' 
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is  described  by  the  vector  13  32  21  12 


21  . 


Thus  the  first  arrow  could  be  kept  fixed  at  the  zeroth 
level  and  the  second  at  the  kth  with  no  difficulty. 


' 
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5 . 6  The  frequency  dependent  susceptibility  series/ 
analysis  and  conclusions 


The  series  we  obtained  for  the  frequency  dependent 
susceptibility  on  the  f.c.c.  lattice  is: 


(u,“2  fi 


I 

s  odd 


s  !/y 

\w 


s+1 


j  X  S+1 


(K) 


X2 (K) =3K+9K2+14 . 75K3+22 . 0K4+61 . 375K5+237 . 8041667K6+ 
880 . 980209K?+  0 (K8) 

x  (K) =20K+6K2-149.5K3-575.291667K4-1433.55K5+  0(K6) 

Xzr  (K)=91.5K+8  3.8  25K2-2  9  6.120  8K3+  0(K4) 

6 

Xg(K)=  391 . 516  66  7K+  0(K2) 


As  part  of  the  calculation  we  had  to  calculate  <[a^(t)  fa^]>/ 
the  bubble  graphs.  As  can  be  seen  from  (  5.5.5)  the  odd 
powers  of  t  for  <[a-^(t)  ,a-^]>  are  just  twice  those  of 
<a^(t)a^>,  so  we  have  the  coefficients  of  the  first,  third, 
fifth  and  seventh  powers  of  t. 

We  also  calculate  the  coefficients  of  the  second  and 
fourth  powers  of  t  in  order  to  see  how  the  dependency  on 

CO 

K  changes  with  the  powers  of  t,  and  estimate  /<a^(t)o^>dt 

0 

(see  (5.2.18)).  The  results  are 


=  1  +  I  a.  (K)  (-i  ait)  i 
i=l 


<a1  (t) a1> 


' 


1 


f 
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a1=3K+6K2+10 .75K3+22.5K4+64„175K5+218.691667K6+ 

+780. 842113 IK 7+  0(K8) 

a2=6+12K+49.5K2+142.0K3+374.25K41054.65K5+3317.85625K6+O(K7) 
a3=28K+97K2+298K3+853.166667K4+2453. 48333K5+0 (K6) 

a  .  =  336  +  1 16 4K+6 70 2K2  +  2 4 39 2K3+0 (K4) 

4 

a  =156 „ 3K+707 . 7K2+2756 . 3125K3+0 (K4) 

5 

a?=664 . 37857K+0 (K2) 


We  assumed  that  <G^(t)a-^>  is  independent  of  N  and 
therefore 


<  q1  (t) g1> 


Tr 


/ 1  \ 

e  G1(t)a1 


Tr 


-3^ 

e 


A  (K  ,  t )  +  NB(K,t)  +  .  .  . 

1  +  a (K)  N  +. . . 


=  A  (K  ,  t) 

A.  Analysis  of  second  and  fourth  moments  of  the  suscep¬ 
tibility 

We  did  the  regular  Pade  approximant  analysis  and  ratio 


analysis  as  seen  in  the  following  tables . 


.  ''S  .1  J 
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TABLE  5.1 


Pade 

approximants 

to  3K  log  *2<K) 

N 

M  =  1 

M  =  2 

M  =  3  M  =  4 

1 

.1754 

.1138 

.2174  .2332 

2 

.1145 

.1479 

.2380 

3 

.2292 

.2547 

4 

.  .2497 

TABLE  5.2 

Pade  approximants  to  (K-Kc)  ^  log  X2 

K= .221 

N 

M  =  0 

M  =  1 

M  =  2 

M  =  3 

M  =  4 

1 

.0407 

.0511 

.0911 

.0925 

.0870 

2 

.0513 

.0377 

.0926 

.0914 

3 

.0987 

.0957 

.0874 

4 

.0954 

.1006 

5 

.0845 

/ 

, 
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TABLE  5.3 
Ratios  of  X2 


n 

rn(X2> 

1 

2 

3 

3 

1.6339 

4 

1.4915 

5 

2.7898 

6 

3. 8746 

7 

3.7046 

TABLE  5.4 

3  _i 

Pade  approximants  to  ^  log  [ (x4 (K) ] 


N  M  =  1  M  =  2 


1  .1891  .2221 


2 


2169 


, 
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TABLE  5.5 


3  “1 

Pade  approximants  to  (K-Kc)  log  X4 (K) 


K= .221 


N 

M  =  0 

M  =  1 

M  =  2 

1 

.735 

.881 

.877 

2 

.859 

r-' 

r- 

00 

• 

3 

.  87  5 

TABLE  5 

.  6 

Pade  approximants  to 

(X4  (K)  ”1) 

1/0 . 8 

N 

M  =  1 

M  =  2 

M  =  3 

M  =  4 

0 

* 

.3481 

.  2443 

.22  35 

1 

* 

.1814 

.2111 

2 

.3199 

.2236 

3 

.2014 

TABLE 

5.7 

Pade 

approximants  (x^ (K) 

1,1/0. 9 

N 

M  =  1 

M  =  2 

M  =  3 

M  =  4 

0 

.3678 

.2561 

.2323 

1 

* 

.1871 

.2167 

3154 


.2282 


3 


2  09  3 
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third 

powers 

of  t  in  <ai(t)ui> 

TABLE 

Ul 

• 

CO 

Pade 

approximants  to  ^ 

log  a1(K) 

N 

M  =  1 

M  =  2 

M  =  3 

M  =  4  M  =  5 

0 

.632 

.414 

.312 

.276  .260 

1 

.352 

.176 

.232 

.240 

2 

.218 

.225 

.242 

3 

.225 

.211 

4 

+ 

.239 

TABLE  5.9 

Pade  approximants  to 

i 

i 

3K  log  al (K)  K=0 .221 

N 

M  =  0 

M  =  1 

M  =  2 

M  =  3  M  =  4  M  =  5 

0 

.4420 

.2679 

.2007 

.1672  .1471  .1333 

1 

.1547 

.0827 

.1007 

.0955  .0852 

2 

.0971 

.0983 

.0968 

.1075 

3 

.0983 

.0976 

.0988 

4 

.0966 

.0989 

.0893 


5 
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TABLE  5.10 


Pade  approximants  to 


a2  (K) 


N 

M  =  1 

M  =  2 

M  =  3 

M  =  4 

M  =  5 

0 

.2131 

.2237 

±.1246i 

.2401 

1 

.4237 

* 

.3333 

.3061 

2 

.7367 

.5149 

.2922 

3 

.2970 

.1277 

4 

.1518 

TABLE  5.11 

Pade  approximants  to  (K-K  )  ^  log  (K) 

K=0 .221 

N  M  =  2  M  =  3  M  =  4 

M  =  5 

3102 

.  3226 

.2475 

.2322 

3218 

.3116 

.2268 

1629 

.0888 

1154 

3 


t 
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TABLE  5.12 


Pade 

approximants 

to  (a2 (K) ) 

1/0.31 

N 

M  =  1 

M  =  2 

M  =  3 

M  =  4 

M  =  5 

1 

.1575 

.1548 

.  1894 
±  .0735i 

.1979 

±.0724i 

.2054 

2 

.1984 

.2736 
±  .0495i 

.2612 
±  ,0601i 

.2803 

3 

.2216 

.2604 
±  .0604i 

.  2695 
±.0545i 

4 

.2389 

.2929 
±  .0346i 

5 

.2507 

TABLE  5.13 

Pade 

approximants 

to  (a2(K)) 

1/.25 

N 

M  =  1 

M  =  2 

M  =  3 

M  =  4 

M  =  5 

1 


2 


3 


4 


140 

* 

.  170± 
.  06  8i 

176 

.224  + 

.  062i 

198 

.  224  + 

.  062i 

.224  + 
.  062i 

216 

.  250± 

.  0 52i 

230 

.192 


5 
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TABLE  5.14 


1/  23 

Pade  approximants  to  ( a£ (K) ) 


N 

M  =  1 

3 

II 

to 

M  =  3 

M  =  4 

M  =  5 

1 

.134 

* 

.  163± 

.  065i 

* 

.186 

2 

.167 

.208  + 

.  062i 

.212  + 

.  061i 

.238 

3 

.190 

.  212± 

.  061i 

.207  + 
.0631 

4 

.207 

.221  ± 

.  053i 

5 

.222 

TABLE  5.15 

Pade  approximants  to  ^  log  a^ (K) 


N 

M  =  1 

M  =  2 

M  =  3 

0 

.373 

.421 

.340 

1 

.415 

.385 

2 


326 


' 


• 
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TABLE  5.16 


Pade 

approximants 

to  (K-K  ) 

w 

d 

3K 

log  a3 (K) 

K=0 .221 

N 

M  =  0 

M  =  1 

M  =  2 

M  =  3 

0 

.766 

.544 

.414 

.341 

1 

.454 

.207 

.139 

2 

.242 

c  118 

3 

.164 

TABLE  5 

.17 

Pade  approximants 

to 

(a3 (K) ) 1/0 

.35 

N 

M  =  1 

M  =  2 

M  =  3 

M  =  4 

101 

.  138± 

.  09  4i 

.221 

.221 

159 

.  197± 

.  0  88i 

.221 

198 

.  244± 

.  075i 

225 

3 


- 

' 
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TABLE  5.18 

1  /  3  8 

Pade  approximants  to  (a^  (K) ) 


N 

M  =  1 

M  =  2 

M  =  3  M  =  4 

0 

.110 

.  155± 

.  lOOi 

.267  .232 

1 

.170 

.213  + 

.001i 

.230 

2 

.209 

.262  + 

.  074i 

3 

.235 

TABLE  5.19 

Ratios 

of  a1 (K) ,  a2 

(K)  and  a^ (K) 

n 

r 

a  i 

r 

a  2 

r 

a3 

1  2 


2 

4.125 

3.4643 

1.7917 

2.8687 

3.0722 

2.0930 

2.6356 

2 . 8630 

2.8522 

2.8180 

2.8757 

3.4077 

3.1459 

3.5705 

7 


- 

! 
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TABLE  5.20 

Pads  approx imants  to  [(K-K^,)0*1  a^  (K)  K  ‘*’]K_  221 


N 

M  =  0 

M  =  1 

M  =  2 

M  =  3 

M  =  4 

M  =  5  M  =  6 

0 

3.920 

3.743 

3.724 

3.726 

3.724  3.715 

1 

3.462 

3.757 

3.721 

3.726 

3.725 

3.728 

2 

3.683 

3.730 

3.727 

3.725 

3.727 

3 

3.722 

3.727 

3.710 

3.728 

4 

3.727 

3.726 

3.729 

5 

3.725 

3.727 

6 

3.720 

- 

TABLE  5.21 

Pade  approximants  to 

[ (K-Kc)0-24  a2(K)]K=221 

N 

M  =  0 

M  =  1 

M  =  2 

M  =  3 

M  =  4 

M  =  5  M  =  6 

0 

5.223 

6.596 

6.495 

5.984 

6.022  6.134 

1 

5.020 

★ 

6.500 

6.619 

6.019 

5.964 

2 

5.879 

6.275 

6.006 

6.132 

6.102 

3 

6.150 

6.114 

6.110 

6.109 

4 

6.109 

6.109 

6.109 

5 

6 . 109 

* 

6 

6.109 

TABLE  5.22 

Pade  approximants  to 

[  (K-K  )°* 

1  c 

35  a3 (K) 

K  ]K=0.221 

N 

M  =  0 

M  = 

=  1 

M  =  2 

M  =  3 

M  =  4 

0 

28 

.25 

26.67 

24.83 

23.77 

1 

23.37 

26 

.78 

34.69 

22.03 

2 

25,65 

25 

.60 

22.31 

3 

25.60 

25 

.65 

4 

24.70 

■ 

The  results  for  the  self  correlation  indicate  that 


<  a 1 (t) a1> 


=  1  + 


n=l 


A2n  (K 


A 

v  \  Y2n  ,  .  ,  ,  2n 

K  )  (-la)  t) 

c  o 


CO 


+  z 

n=0 


A2n+1 


K  (K 


A 

^2n+l  ,  .  ,  v 2n+l 

(no  t) 


with 


y^  =  0.099  ±  0.005 
yA  =  0.24  ±0.03 

yA  =  0.35  ±  0.03 


A±  =  0.82  ±  0.1 
A^  =  6.1  ±  0.5 
A3  =  5.3  ±2 


(5.6.1) 


where  the  A  were  obtained  by  taking  Pade  approximants 
^  yA 

of  the  function  [(K-K  )  n  a^  (k|  evaluated  at  K=Kc  for  n 

YA  —  1 

and  of  [ (K-K  )  n  a  (K)K  ]  evaluated  at  K=K  for  n  odd. 
on  ^ 

A  A 

The  differences  are  A  =  y  -  y  -  0.10  ±  0.05 

xr  11  II  X 


-n  A 

<a1(t)a1>  -  1  +  Z  An(T-Tc)  Afiwot)n 


'  t  ' 


A 


(5.6.2) 


Near  T  if  we  can  assume  that  f  f(x)dx  converges 
c'  0 


A 


we 


A 

find  that  the  autocorrelation  will  vanish  like  e 


<  a,  (t)  a,  >  ~  /  f 

l  i  o 


A 


A 


dt  =  £  A/f(x)dx  -  e  A 


(5.6.3) 


CO 

We  saw  that  A^MR  /  <0-^  (t )  a^>dt  -  e 


(5.6.4) 


and  in  (5.2.18)  we  obtained  that  A 

2 

in  our  case  n  =  0 ,  v  =  ±  , 


NMR 


-  £ 


v  ( n+i-z) 


even 
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therefore  z  =  0.85  ±  0.10.  As  seen  in  (5.2.20) ,  z  -  ; 

therefore  A  =  0.57  ±  0.10 ,  where  A  is  the  critical  slow- 
s  s 

ing  down  index. 

The  result  for  A,  is  corroborated  by  the  calculation 

ri 

of  A  from  the  frequency  dependent  susceptibility  and  then 
s 

substituting  in  (5.2.18)  to  obtain  A^. 

As  we  know 

— v 

X  (t0=0  ,£)-£. 

The  behaviour  in  oo  we  shall  assume  to  be 

-y/A 

t  £=0)  -  a)  S  ,  (5.6.5) 


where  this  defines  A  . 

s 

A  possible  function  with  this  behaviour  and  only 
even  moments  is 


X  (w) 


2A 

s 

+  £ 


f 


-Y/2As 


(5.6.6) 


and  f  ( oo  /  £ )  regular  at  oo=0  ,  e=0  . 

The  index  of  critical  slowing  down  defined  by  Suzuki 
(1969)  as  mentioned  in  (5.1.10)  is 


i.  = 


lim 

co->0 


ym  -  x(or 


=  lim 
cu-v  0 


2 

) 

(1-f 


00  x(0) 

-1 


,2  -Y/2A 

to  .  s 


-  1 


00 


-A 


=  £ 


(5.6.7) 


i 
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Our  assumption  is  that  x^)  scales  like  f 


e 


therefore  the  moments 


“Y 4  -fy-4A  > 

X  ~  e  -  e 

4 


From  our  results 

y  =  +.096  ±  .005  -*  A  =  .62  ±  .10 
2  s 

y  =  -.88  ±  .05  +  A  =  .55  ±  .10  . 

1  4  s 

This  agrees  with  our  former  estimate. 

Clearly  our  estimates  of  A  however  large  the  error 

s 

margins  yield 


since 


0.58  ±  0.10  . 


A 


s 


This  agrees  with  Tomita  (1968)  for  the  isotropic  Heisenberg 
model  where  kinetic  speeding  up  occurs,  while  the  kinetic 
Ising  model  had  A  >  y  a  kinetic  slowing  down. 

S  ' 


■ 

t 

■ 


' 


* 


A1 


APPENDIX  A 


Graphs  and  Weights 

In  the  following  tables  the  graphs  needed  for  the 
calculation  of  the  static  initial  susceptibility, 
dynamical  susceptibility  and  fourth  order  fluctuation 
are  listed. 

Table  A .5  lists  the  fluctuation-like  graphs  used 
for  the  static  susceptibility.  The  bubble  graphs  are 
not  listed  as  they  are  identical  to  those  of  the  dynamical 
susceptibility  listed  in  Table  A. 7  apart  from  the  vertical 
weight.  The  static  weight  is  the  sum  of  all  dynamical 
weights.  Table  A. 9  has  the  partition  function  graphs. 

Table  A. 6  lists  the  fluctuation -like  graphs  generated 
for  the  dynamical  susceptibility. 

Table  A. 8  lists  the  fourth  order  fluctuation  graphs. 

The  entries  are  in  two  lines  for  each  graph.  In 
the  first  line  in  order  the  entries  are:  the  number  of 
vertices,  the  number  of  solid  arrows  (except  for  Table 
A. 2  which  lists  the  number  of  all  arrows,  solid  plus  two 
dotted) ,  lattice  constant,  horizontal  weight,  and  vertical 
weight  or  weights.  The  lattice  constants  are  the  f.c.c.  ones. 

In  the  second  line  the  graph  itself  is  listed,  for 
example 

1  223311221  . 


, 


A  2 


Dotted  arrows  come  first  in  the  listing.  In  Tables  A . 6 , 
A. 7  the  spurious  vertex  is  kept  in  the  listing  of  the 
graph,  for  example 


2442233112  . 

Since  v(k)  =  v(£-k)  we  did  not  list  all  vertical 
weights,  only  4  for  £  ^  7  and  5  for  £  =  8  . 

Below  is  the  subdivision  of  graphs  according  to 
the  number  of  solid  arrow  £. 

TABLE  A . 1 

Fluctuation-like  static  susceptibility  graphs 


£ 

No.  of  graphs 

1 

1 

2 

1 

3 

4 

4 

9 

5 

30 

6 

80 

7 

268 

8 

854 

Total 

1247 

A3 


TABLE  A. 2  TABLE  A. 3 

Fluctuation-like  dynamic  Bubble  graphs 

susceptibility  graphs 


£ 

No.  of  graphs 

£ 

No.  of  graphs 

3 

1 

2 

1 

4 

2 

3 

1 

5 

10 

4 

5 

6 

26 

5 

8 

7 

103 

6 

32 

8 

348 

7 

69 

8 

245 

Total 

4  90 

Total 

361 

TABLE  A . 4 

Fourth  order  fluctuation  graphs 


£ 

No.  of  graphs 

2 

2 

3 

3 

4 

20 

5 

61 

6 

292 

7 

1117 

Total 

1495 
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